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Vitreous  silica  is  a material  of  considerable  techni- 
cal and  scientific  importance  and  hence  its  structure  and 
properties  have  been  widely  studied.  The  two  most  important 
contributions  to  the  study  of  the  structure  of  vitreous 
silica  were  made  by  Warren  working  with  different  co-workers 
at  two  different  periods  separated  by  about  30  years.  In 
1936,  Warren,  Krutter  and  Morningstar  published  an  approxi- 
mate analysis  of  the  structure  of  vitreous  silica.  They 
were  able  to  obtain  the  Radial  Distribution  Function  of 
vitreous  silica  and  show  that  their  results  were  consistent 
with  the  random  network  theory  of  the  structure  of  vitreous 
silica.  In  1969,  Mozzi  and  Warren  published  the  results  of 
a new  structural  study  of  vitreous  silica.  Using  the 
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precise  Pair  Function  Distribution  (PFD)  approach,  they  were 
able  to  analyze  equally  precise  x-ray  intensity  data  to  ob- 
tain a better  description  of  the  local  atomic  distributions 
in  vitreous  silica.  They  used  Distance  Distribution  Func- 
tions (DDF)  to  represent  the  distribution  of  atomic  dis- 
tances in  vitreous  silica. 

This  dissertation  presents  and  applies  an  approach  for 
using  the  DDF  to  characterize  structural  changes  occurring 
in  vitreous  silica  as  a function  of  heat- treatment , water 
content, and  as  a result  of  fast  neutron  irradiation.  Bulk 
and  powder  samples  of  "T-D"  clear  fused  quartz  were  heat- 
treated  at  different  temperatures  and  the  resulting  struc- 
tural changes  studied  by  an  x-ray  diffraction  analysis.  It 
was  found  that  the  heat -treatment  of  vitreous  silica  does 
not  produce  any  significant  changes  in  the  bulk  structure 
of  vitreous  silica.  The  powder  heat-treatment  results  sug- 
gest that  the  structure  of  vitreous  silica  may  be  undergoing 
some  changes  during’  the  initial  stages  of  surface  crystal- 
lization. The  DDF  analysis  has  been  applied  to  follow  the 
structural  changes  in  powdered  vitreous  silica  as  a func- 
tion of  time  at  heat -treatment  temperature . The  effect  of 
water  content  on  the  structure  of  vitreous  silica  was  stud- 
ied using  samples  of  IR  Vitreosil,  Vitreosil  066  and 
Spectrosil-B.  These  samples  span  a range  of  water  contents 
from  0.0003  percent  in  the  case  of  IR  Vitreosil  to  0.12  per- 
cent. in  the  case  of  Spectrosil-B.  It  was  found  that  the 
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PFD  of  IR  Yitreosil  and  Vitreosil  066  correspond  closely  to 
the  PFD  of  "T-D"  clear  fused  quartz,  while  the  PFD  of  Spec- 
trosil-B  indicates  that  the  structure  of  the  latter  may  be 
more  "crystalline"  than  that  of  the  other  samples.  The 
present  study  also  shows  that  the  PFD  of  an  Optosil  sample 
undergoes  a significant  change  as  a result  of  heavy  irradi- 
ation by  fast  C > 1 0 KeV)  neutrons.  The  DDF  analysis  of  this 
sample  reveals  considerable  changes,  especially  in  the  near- 
est neighbor  distributions  as  a result  of  neutron  irradia- 
tion. 
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CHAPTER  I 


INTRODUCTION 

•Vitreous  silica  is  a material  of  considerable  impor- 
tance' from  a technological  as  well  as  scientific  point  of 
view..  It  is  the  only  commercial  one-component  glass,  and 
has  many  applications  owing  to  its  high  softening  point,  re- 
sistance  to  chemical  attack  and  its  high  transparency  in  the 
visible  and  ultra-violet  regions.  Scientifically  it  is  im- 
portant as  being  the  simplest  member  of  the  family  of  sili- 
cate glasses.  Hence  it  is  not  surprising  that  the  struc- 
ture and  properties  of  vitreous  silica  have  been  widely 
studied.  The  systematic  study  of  the  structure  of  vitreous 
silica  started  as  early  as  the  1930's.  Warren,  Kruttcr 
and  Morningstar  (1936)  were  the  first  to  analyze  the  x-ray 
diffraction  pattern  of  vitreous  silica  to  obtain  its  Radial 
'Distribution  Function  (RDF).  However,  their  analysis  was 
only  approximate  as  they  had  to  assume  that  the  scattering 
•factors  of  silicon  and  oxygen  were  proportional  to  one 
another  in  order  to  successfully  carry  out  the  Fourier  analy- 
sis of  the  vitreous  silica  diffraction  data.  This  assump- 
tion is  not  strictly  correct.  Also,  no  allowance  was  made 
to  account  for  the  termination  ripples  produced  by  the 
Fourier  analysis  of  the  diffraction  data. 
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In  the  1950's,  a new  theoretical  approach  was  devel- 
oped to  treat  the  diffraction  data  from  glasses  such  as 
vitreous  silica  without  the  need  for  any  approximation. 

This  was  the  pair  function  analysis.  The  pair  function 
analysis  takes  into  account  the  fact  that  in  a glass  con- 
taining more  than  a single  element  the  atomic  scattering 
factors  of  the  constituent  atoms  need  not  be  proportional 
to  one  another.  The  pair  function  analysis  was  applied  to 
the  case  of  vitreous  silica  by  Mozzi  and  Warren  (1969) . 

They  were  able  to  measure  accurate  intensity  data  to  high 
values  of  S(2  s^-n-— ) and  obtain  a well -resolved  Pair 
Function  Distribution  (PFD)  of  vitreous  silica.  The  inter- 
pretation of  the  PFD  takes  the  termination  ripples  into  ac- 
count. Mozzi  and  Warren  described  their  results  in  terms 
of  Distance  Distribution  Functions  (DDF)  which  are  related 
to  the  statistical  distribution  of  interatomic  distances  in 
vitreous  silica.  The  PFD-DDF  approach  gives  a precise  quan- 
titative statistical  description  of  the  structure  of  vitre- 
ous silica. 

In  the  present  investigation  the  PFD-DDF  approach  has 
been  used  t.o  study  the  structural  changes  in  vitreous  silica 
as  a function  of  heat-treatment , water  content  and  heavy 
neutron  irradiation. 


CHAPTER  II 


STRUCTURE  OF  VITREOUS  SILICA 

Just  as  sound  is  characterized  by  its  Fourier  spectrum, 
so  a state  of  aggregation  of  a substance  is  most  simply 
characterized  by  the  Fourier  transform  of  its  scattering 
power,  i.e.,  by  the  Fourier  transform  of  its  x-ray,  neutron 
or  electron  diffraction  pattern  [Prins  (1968)].  The  Fourier 
transforms  of  diffraction  patterns  can  yield  the  RDF  (or  the 
PFD)  of  the  glass.  Consequently,  the  diffraction  patterns 
or  the  RDF  (PFD)  can  be  used  to  characterize  a glass  like 
vitreous  silica.  The  RDF  gives  a description  of  the  struc- 
ture of  glass,  i.e.,  the  distribution  of  atoms  in  the  glass. 

Survey  of  Pre vious  Work  on  the 
Structure  of  Vitreous  Silica 

The  first  systematic  investigation  of  the  structure 
of  vitreous  silica  was  the  x-ray  study  of  vitreous  silica 
by  Randall,  Rooksby  and  Cooper  (1930).  They  found  that  the 
x-ray  diffraction  patterns  of  vitreous  silica  contain  broad 
maxima  at  the  same  positions  as  the  strongest  lines  in  crys- 
talline silica.  They  concluded  that  vitreous  silica  was  an 
aggregate  of  extremely  small  crystals  (crystallites).  A 
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completely  different  view  of  the  structure  of  vitreous 
silica  was  proposed  by  Zachariasen  (1932)  . Zachariasen  sug- 
gested that  vitreous  silica  consisted  of  a three-dimensional 
network  of  randomly  oriented  SiO^  tetrahedra. 

Random  Network  Theory  vs.  the  Crystallite  Hypothesis 

The  first  radial  distribution  analysis  of  the  structure 
of  vitreous  silica  was  that  of  Warren,  Krutter  and  Morning- 
star  (1S36)  . Using  powder  cameras  to  record  the  diffracted, 
intensities,  they  arrived  at  an  RDF  of  vitreous  silica. 

From  their  experimental  RDF  they  concluded  that  the  struc- 
ture of  vitreous  silica  could  indeed  be  described  as  a con- 
tinuous random  network  of  Si.O^  tetrahedra.  Warren  (1937) 
also  conducted  x-ray  small  angle  scattering  studies  of  vit- 
reous silica.  He  could  not  detect  any  small  angle  scatter- 
ing from  vitreous  silica  and  hence  concluded  that  vitreous 
silica  could  not  be  an  aggregate  of  small  crystals. 

Valenkov  and  Porai - Koshits  published  the  results  of 
their  x-ray  studies  of  vitreous  silica  in  1936. . They  sug- 
gested that  vitreous  silica  need  not  consist  of  a continuous 
random  network.  Individual  molecules  or  ions  in  the  glass 
could  be  assembled  into  large  or  small  groups.  The  order 
within  the  groups  \vill  be  larger  than  the  mean  order  in  the 
glass  and  the  groups  can  be  regarded  as  small  distorted 
crystals  separated  from  one  another  by  zones  of  lesser  order. 
This  is  the  so-called  crystallite  hypothesis  of  the  struc- 
ture of  glass. 
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Hartlief  s Analysis  of  the  Structure  of 
V it reous  Silica 

Warren  and  Porai-Koshits  have  both  used  x-ray  tech- 
niques to  arrive  at  different  conclusions  about  the  struc- 
ture of  vitreous  silica.  In  an  attempt  to  resolve  these 
differences,  Hartlief  (1958)  undertook  a careful  x-ray  study 
of  the  structure  of  vitreous  silica.  Hartlief  used  ioniza- 
tion chambers  to  record  the  diffracted  intensity.  He  ob- 
tained an  RDF  which  agreed  well  with  the  RDF  obtained  by 
Warren  and  co-workers  except  for  the  presence  of  two  addi- 
tional peaks  at  6.3  A and  7.5  A.  Hartlief  suggested  that 
the  difference  between  his  and  Warren's  RDF  might  be  due  to 
the  differences  in  the  experimental  techniques,  or  due  to 
inherent  differences  in  the  structure  of  the  samples  stud- 
ied. Hartlief  can  be  said  to  be  the  first  investigator  to 
consider  the  possibility  that  the  structure  of  vitreous 
silica  might  be  influenced  by  its  history. 

S i - 0-S i Bond  Angle  in  Vitreous  Silica 

The  description  of  the  distribution  of  the  Si-O-Si 
bond  angle  in  vitreous  silica  is  of  importance  in  character- 
izing the  structure  of  vitreous  silica.  Zarzycki  (1957) 
has  investigated  this  question  rather  thoroughly.  From 
structural  and  energetic  considerations,  Zarzycki  concluded 
that  the  minimum  possible  Si-O-Si  angle  in  vitreous  silica 
is  116°,  and  that  the  angle  varies  from  116°  to  180°  with 
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a median  value  of  about  148°.  Zarzycki  compared  this  dis- 
tribution to  the  Si-O-Si  bond  angle  in  low-cristobalite , 
which  is  150°. 

Other  Investigations  of  the  Structure 
bT~~V:i  treous  Silica 

Simon  (1960)  and  more  recently  Bale,  Shepler  and 
Gibbs  (1970)  have  also  used  the  Warren  analysis  to  obtain 
the  RDF  of  vitreous  silica. 

The  investigations  we  have  considered  so  far  were  all 
based  on  x-ray  diffraction  studies.  Neutron  diffraction 
methods  have  also  been  used  to  study  the  RDF  of  vitreous 
silica  [Breen,  Delaney,  Persiani  and  Weber  (1957);  Delaney 
and  Weber  (i960);  Carraro,  Domenici  and  Zucca  (1965);  and 
Lorch  (1969) ] . 

Keating  (1963)  suggested  that  a combination  of  neutron 
and  x-ray  diffraction  methods  can  be  used  to  interpret  the 
structural  details  of  a liquid-like  glass.  Henninger, 
Buschert  and  Heaton  (1967)  obtained  vitreous  silica  RDF 
using  both  x-ray  and  neutron  diffraction  data.  The  theo- 
retical analysis  was  again  the  same  as  that  used  by  Warren, 
Krutter  and  Morningstar  (1936).  The  PFD  method  was  first 
applied  by  Mozzi  and  Warren  (1969)  to  study  the  structure 
of  vitreous  silica.  The  PFD  analysis  is  more  precise  than 
the  RDF  method.  Mozzi  and  Warren  were  able  to  obtain  a 
well  - resolved  picture  of  the  structure  of  vitreous  silica. 
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They  were  able  to  give  a precise  quantitative  description 
of  the  variation  of  the  Si-O-Si  bond  angle  as  well  as  the 
variation  of  the  orientation  of  SiO^,  tetrahedra  about  Si-0 
bond  directions.  They  concluded  that  the  results  of  the 
PFD  analysis  were  consistent  with  a random  network  model  of 
vitreous  silica. 

Fine  Structure  in  Vitreous  Silica 

It  is  also  of  interest  to  know  how  the  SiO^  tetrahedra 
are  joined  together,  i.e.,  to  know  how  many  of  the  rings  in 
the  SiO^  network  are  five,  six  or  seven  membered.  Oberlies 
and  Dietzel  (1957)  and  Robinson  (1965)  have  attempted  to 
postulate  models  to  describe  the  distribution  of  ring  sizes 
in  vitreous  silica.  The  diffraction  methods  cannot  be  ex- 
pected to  give  the  details  of  this  type.  What  the  diffrac- 
tion methods  can  reveal  is  a picture  which  is  a statistical 
average  of  the  local  atomic  arrangements  in  vitreous  silica. 

Construction  of  Models 

A statistical  description  of  the  structure  of  vitreous 
silica  can  also  be  obtained  by  construction  of  ball  models. 
Ball  and  wire  models  are  wridely  used  in  the  study  of  crys- 
tal structures.  Ordway  (1964),  Evans  and  King  (1966),  and 
Bell  and  Dean  (1966)  have  applied  this  technique  to  study 
the  structure  of  vitreous  silica.  These  investigators 
used  essentially  the  following  procedure.  A model  was 
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constructed  by  joining  together  SiO^  tetrahedron  units  in 
an  essentially  random  fashion.  The  distribution  of  bond 
lengths  as  \\'ell  as  the  average  number  of  neighbors  could 
be  measured  from  the  model.  Their  results  are  essentially 
similar  to  the  RDF. 

Summary  of  Previous  Wor k on  the 
Structure  of  Vitreous  Silica 

The  conclusions  of  the  various  investigators  as  to  the 
structure  of  vitreous  silica  are  summarized  in  Table  1.  It 
should  be  noted  that  most  of  the  investigators  have  reported 
mainly  the  interatomic  distances  in  vitreous  silica. 

Zarzycki  has  given  the  possible  variation  in  the  Si-O-Si 
bond  angles,  while  Mozzi  and  Warren  have  reported  a func- 
tion V(a)  which  gives  a precise  quantitative  description  of 
this  variation. 


Influence  of  History  on  the 
Structure  ~o~f  Vitreous  Silica 

As  we  have  seen,  the  structure  of  vitreous  silica  has 
been  widely  studied.  The  same  is  true  of  the  properties 
of  vitreous  silica  as  well.  It  is  now  well  established 
that  the  properties  of  vitreous  silica  are  affected  by  vari- 
ables such  as  thermal  history,  impurity  content,  and  treat- 
ments such  as  neutron  irradiation.  However,  except  for  the 


Interionic  Distances  in  Vitreous  Silica  (A) 
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case  of  neutron  irradiation,  there  have  been  no  systematic, 
structural  studies  conducted  as  a function  of  these  vari- 
ables. 

Ef fee t of  He at -Treatment 

For  example,  many  studies  have  been  made  of  the  crys- 
tallization kinetics  in  vitreous  silica  [Ainslie  et  al . 
(1962),  Wagstaff  et  al . (1964)].  One  of  the  more  interest- 
ing studies  from  the  structural  point  of  view  is  that  of 
Wagstaff  (1968)  on  the  growth  rate  of  internally  nucleated 
cristobalite  in  vitreous  silica.  The  growth  of  cristobalite 
in  bulk  glass  can  be  expected  to  alter  the  structure  of  the 
surrounding  glass.  It  would  seem  that  by  choosing  proper 
experimental  conditions  it  would  be  possible  to  study 
structural  changes,  if  any,  in  the  bulk  of  vitreous  silica 
as  it  becomes  crystalline.  The  PFD  method,  for  example, 
can  be  applied  to  study  such  a structural  change.  This  is 
the  basis  of  part  of  the  present  investigation. 

Zarzycki  (1957)  studied  the  x-ray  diffraction  patterns 
of  vitreous  silica  at  20°  and  1600°C.  He  obtained  the  RDF 
at  both  the  temperatures  and  found  essentially  no  change 
in  the  structure  of  vitreous  silica  at  the  two  temperatures 
studied.  However,  Zarzyc'k.i  did  not  determine  if  any  changes 
occurred  in  the  structure  of  vitreous  silica  as  a function 
of  time  at  the  elevated  temperature. 
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The  crystallization  kinetics  cf  powdered  vitreous 
silica  have  been  studied  using  standard  x-ray  diffraction 
procedures.  Verduch  (1958)  and  Coquerelle  (1961)  con- 
cluded that  an  intermediate  stage  could  be  detected  before 
the  powdered  vitreous  silica  became  crystalline  at  temper- 
atures below  1300 °C,  while  March and  and  Favede  (1968)  could 
detect  no  such  intermediate  stage  in  their  investigations 
of  the  crystallization  kinetics  of  a sample  of  powdered 
vitreous  silica. 

Effect  of  Water  Content 

There  has  been  no  systematic  study  of  the  effect  of 
water  content  on  the  structure  of  silica.  This  is  surpris- 
ing because  the  effect  of  water  content  on  properties  of 
vitreous  silica  such  as  viscosity  [Hetherington,  Jack  and 
Kennedy  (1964)],  density  [Bruckner  (1970)],  and  crystalli- 
zation kinetics  [Wagstaff  et  al . (1964)]  has  been  investi- 
gated and  found  to  be  considerable. 

Effect  o f Neutron  Irradiation 

On  the  other  hand,  RDF  methods  have  been  used  to  in- 
vestigate the  effect  of  neutron  irradiation  on  the  struc- 
ture of  vitreous  silica.  Primak  and  co-workers  were  among 
the  first  to  study  the  effect  of  neutron  irradiation  on  the 
properties  of  vitreous  silica  [Primak  (1958),  Primak  and 
Kampwirth  (1968)].  Stolen,  Krause  and  Kurkijian  (1971) 
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have  observed  significant  changes  in  the  optical  and  thermal 
properties  of  vitreous  silica  as  a result,  of  heavy  neutron 
irradiation.  Simon  (1960)  and  Lorch  (1969)  ha\?e  studied 
the  change  in  the  structure  of  vitreous  silica  as  a result 
of  heavy  neutron  irradiation.  Bale,  Shepler  and  Gibbs  (1970) 
have  studied  the  structural  changes  in  a GE  vacuum  fused 
vitreous  silica  as  a function  of  various  dosages  of  neutron 
irradiation.  The  results  of  these  investigations  can  be 
summarized  as  follows. 

1.  The  structure  of  heavily  irradiated  vitreous  silica 
is  identical  to  the  structure  of  heavily  irradiated  crystal- 
line quartz  [Simon  (I960)]. 

2.  The  mean  value  of  the  Si-O-Si  bond  angle  shifts 
to  lower  angles  with  increasing  neutron  irradiation. 

3.  The  RDF  peaks  of  vitreous  silica  become  broader 
with  increasing  irradiation,  indicating  a general  displace- 
ment of  atoms  from  their  original  positions. 

The  present  work  is  a study  of  the  structural  changes 
in  vitreous  silica  as  a result  of  heat -treatments , changes 
in  water  content  and  heavy  neutron  irradiation. 


CHAPTER  III 


X-RAY  DIFFRACTION  ANALYSIS 
OF  VITREOUS  STRUCTURES 

Introduction 

Debye  (1915)  and  Ehrenfest  (1915)  independently  sug- 
gested that  crystallinity  is  not  an  essential  prerequisite 
for  the  production  of  an  x-ray  diffraction  pattern.  They 

also  suggested  that  the  x-ray  diffraction  pattern  of  a non- 

* 

crystalline  substance  should  be  characteristic  of  the  sub- 
stance and  can  be  used  to  study  the  atomic  distributions  in 
the  non- crystalline  substance.  Zernike  and  Prins  (1927) 
proposed  that  the  x-ray  diffraction  data  from  monoatomic 
liquids  can  be  treated  by  Fourier  transform  methods  to 
yield  the  average  atomic  distributions  in  these  liquids. 
Debye  and  Menke  (1930)  were  the  first  to  apply  Fourier 
transform  methods  to  study  the  structure  of  a liquid.  .Sub- 
sequently, Warren  and  Gingrich  (1934)  used  a modified  form 
of  this  analysis  to  study  the  structure  of  various  liquids 
and  glasses.  This  chapter  contains  a discussion  of  the  ap- 
plication of  x-ray  diffraction  methods  to  the  analysis  of 
the  structure  of  non-crystalline  solids  with  special  refer- 
ence to  vitreous  silica. 
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X-ray  Diffraction  by  Non- crystalline  Solids 

Let  us  consider  a collimated  beam  of  monochromatic 
x-rays  of  wavelength  y being  scattered  by  a piece  of  glass 
(Fig.  1).  The  atoms  which  constitute  the  glass  will  scatter 
the  x-rays  in  all  directions.  The  scattered  x-rays  along 
any  particular  direction  making  an  angle  of  28  with  the  in- 
cident beam  would,  in  general,  consist  of  two  components: 
one  the  unmodified,  or  coherent,  and  two  the  modified,  or 
the  incoherent,  component.  The  modified  component  does  not 
contain  any  direct  structural  information  and  will  be  ig- 
nored in  the  following  discussion. 

The  average  unmodified  scattering  along  any  direction 
29  is  given  by 


I(S) 


sin  27rSr 


E Z f f 
m „ m n 
m n 


nm 


2irSr 


nm 


(3.1) 


where  I (S)  is  the  unmodified  scattered  intensity  from  the 
sample  in  absolute  units  neglecting  absorption,  f and  f 
are  the  scattering  factors  of  atoms  "m"  and  "n"  respectively, 
"m"  and  "n"  denote  the  atoms  separated  by  a distance  of 
r A,  and  "m"  and  "n,!  can  refer  to  the  sarnie  or  different 
atomic  species. 

Equation  (3.1)  was  originally  developed  by  Debye  (1913) 
and  is  generally  referred  to  as  the  Debye  equation.  An 
examination  of  the  Debye  equation  reveals  that  the  unmodi- 
fied scattering  from  a glass  contains  information  regarding 
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Figure  1.  X-ray  diffraction  from  a non-crystalline 
sample  (schematic) . 
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the  interatomic  distances,  r , the  average  number  of  atoms 
at  distances  r , and  the  type  of  atoms  at  these  distances. 
In  other  words,  the  Debye  equation  contains  information  re- 
garding the  structure  of  glass.  Any  study  of  the  structure 
of  a glass  can  be  then  divided  into  two  important  steps. 

The  first  step  is  to  obtain  the  unmodified  scattering  in- 
tensity data  from  the  glass  as  a function  of  S,  and  the 
second  step  is  to  analyze  the  unmodified  intensity  data  so 
as  to  obtain  the  structural  information.  The  structural  in- 
formation obtained  in  this  manner  is  often  represented  in 
the  form  of  an  RDF. 

Radial  D istribution  Functions 

For  a monoatomic  glass,  the  Debye  equation  can  be  di- 
rectly used  to  treat  the  intensity  data  to  give  the  RDF  of 
the  glass.  The  RDF  of  a glass  gives  the  atomic  density 
distribution  about  any  atom  in  the  glass.  The  RDF  is  gen- 
erally shoim  plotted  as  a function  of  the  distance  from  a 
central  atom  and  consists  of  peaks  and  dips  corresponding 
to  the  fluctuations  in  the  atomic  density  distributions 
about  the  central  atom.  The  peak  positions  represent  the 
mean  interatomic  distances  in  the  glass  and  the  area  under 
each  peak  corresponds  to  the  number  of  atoms  situated  about 
that  mean  distance  of  separation.  The  peaks  and  dips  die 
down  after  a distance  of  about  10  X.  When  the  glass  under 
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study  contains  more  than  one  atomic  species  the  Fourier 
transform  of  the  diffraction  data  is  not  straightforward. 

The  Radial  Pi s t ribution  Analysis 
oT~Vitfe6us  Silica 

Vitreous  silica  consists  of  two  different  atomic  species, 
i.e.,  silicon  and  oxygen.  Hence  the  vitreous  silica  dif- 
fracted intensity  data  cannot  be  analyzed  in  a straightfor- 
ward manner.  Warren,  Krutter  and  Mornings tar  (1936)  were 
the  first  to  suggest  a solution  to  this  problem.  They  as- 
sumed that  the  scattering  factors  of  the  silicon  and  oxygen 
atoms  are  proportional  to  one  another.  By  virtue  of  this 
assumption,  the  vitreous  silica  intensity  data  become  amen- 
able to  Fourier  transform  methods  and  the  vitreous  silica 
RDF  could  be  obtained.  This  method  was  only  approximate  be- 
cause of  the  following  factors: 

1.  The  assumption  of  proportionality  between  the  atom- 
ic scattering  factors  of  silicon  and  oxygen  is  not  strictly 
correct  [Warren  (1969a)]. 

2.  Termination  ripples  which  occur  during  the  Fourier 
transform  of  the  intensity  data.  X-ray  diffraction  data  are 
only  available  up  to  a finite  value  of  S,  thus  limiting  the 
upper  limit  of  the  Fourier  transform  integral.  The  limiting 
of  the  Fourier  transform  integral  to  a finite  S value  intro- 
duces ripples  throughout  the  RDF.  These  ripples  can  cause 
considerable  confusion  in  the  interpretation  of  the  RDF. 
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These  two  shortcomings 
by  the  use  of  the  Pair 


of  the  RDF  method  can  be  eliminated 
Function  Distribution  (PFD)  method. 


Pair  Function  Distribu t i on 

Waser  and  Schomaker  (1953)  pointed  out  that  rigorous 
methods  exist  for  the  analysis  of  diffraction  data  from 
glass  samples  containing  more  than  one  kind  of  atom.  Finbak 
and  co-workers  have  developed  the  rigorous  methods  and  their 
treatment  yields  [Warren  (1969a,b)]„ 

N.  . 
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where  u.c.  denotes  one  unit  of  composition,  SiC>2 , for 
vitreous  silica,  and  k = 2 tt S . 
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= The  unmodified  scattering  intensity  at  large 
angles  in  electron  units  per  unit  of  compo- 
sition. 

= A sharpening  factor  which  is  chosen  to  be 
approximately  unity  at  k = 0 . 
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The  atomic  scattering  factor  for  atom  "j" 
corrected  for  anomalous  dispersion,  a function 
of  k. 

The  number  of  electrons  in  atom  "jM. 

The  number  of  neighbors  in  the  ith  shell  about 
atom  ’'j”. 

The  distance  from  an  atom  "j"  to  an  atom  in 
the  i^*1  shell. 

The  cut-off  value  of  k. 

A convergence  factor  introduced  to  minimize 

termination  satellites  and  weight  down  the 

less  accurate  measurements  at  high  k values, 

ak  ~ 1 . 
m 

The  convolution  of  a bond  length  distribution 
function  with  the  pair  function  P.  . (r ) given  by 
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P.j(r)  can  be  easily  calculated  by  using  the  auxiliary  pair 
functions : 
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and  the  equality 
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P-.(r)  ~ Q. -(r-r..)  + Q..(r+r..) 
a 3 13  ij  13  13 


(3.6) 


with  Qj[j(r+rjj)  usually  negligible. 

Equation  (3.2)  gives  the  PFD  of  the  atoms  in  the  sam- 
ple under  study.  The  experimental  intensity  values  are 
used  to  calculate  the  Fourier  integral  in  the  right-hand 
side  of  Eq . (3.2).  The  right-hand  side  of  Eq.  (3.2)  gives 
the  experimental  PFD.  The  left-hand  side  of  Eq . (3.2)  is 
then  calculated  assuming  various  structural  models  for  the 
glass  under  study  to  give  calculated  PFD  until  a good  fit 
is  obtained  between  the  calculated  and  experimental  PFD  for 
a particular  structural  model.  This  procedure  makes  it  pos- 
sible to  describe  the  structure  of  glass  in  a more  precise 
manner  than  was  possible  by  the  old  RDF  method.  Mozzi  and 
Warren  (1S69)  applied  the  PFD  method  to  analyze  the  struc- 
ture of  vitreous  silica.  Using  sophisticated  experimental 
techniques  they  were  able  to  obtain  accurate  intensity  data 
up  to  k = 20  A They  analyzed  the  intensity  data  using 

the  PFD  method  and  their  results  are  shown  in  Fig.  2.  The 
solid  curve  in  Fig.  2 is  the  experimental  PFD  obtained  by 
Mozzi  and  Warren  for  vitreous  silica.  The  calculated  PFD 
is  shown  as  a dashed  line.  The  theoretical  approach  of 
Mozzi  and  Warren  can  be  summarized  as  follows. 
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Figure  c. . Comparison  of  experimental  and  calculated  PFD  for  vitreous 
silica  [adapted  from  Mozzi  and  Warren  (1969)]. 
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Pair  Function  Analysis  of  the 

— ■ — - - - - * — 

Structure  of  Vitreous  Silica 

In  crystalline  silicates  it  is  known  that  each  silicon 
is  tetrahedrally  surrounded  by  four  oxygens  with  a silicon- 
oxygen  distance  of  about  1.62  A and  an  oxygen- oxygen  dis- 
tance of  about  . 2. 65  A.  With  this  information  the  peaks  in 
the  vitreous  silica  PFD  (Fig.  2)  can  be  identified  as  fol- 
lows. The  first  peak  in  the  PFD  is  due  to  the  interaction 
of  all  the  silicon  and  oxygen  atoms  at  a distance  of  separa- 
tion of  about  1.62  A per  unit  of  composition.  Ths  will  be 
called  a Si-1  0 interaction.  Similarly,  the  second  peak 

at  2.65  A is  due  to  an  oxygen-lst  oxygen  (0-lst  0)  inter- 

s t s t 

action,  the  third  peak  due  to  silicon- 1 silicon  (Si-1  Si) 
interaction,  the  fourth  due  to  a silicon- 2n<^  oxygen  (Si-2nc* 

0)  , and  the  fifth  peak  due  to  a combination  of  oxygen-2IU 
oxygen  (0-2n<^  0)  and  silicon- 2nc^  silicon  (Si-2UC^  Si)  inter- 
actions per  unit  of  composition. 

Having  identified  the  peaks  of  the  PFD,  the  next  step 
is  to  determine  the  distribution  of  the  interatomic  distances 
in  vitreous  silica.  This  distribution  of  atomic  distances 
can  theoretically  be  calculated  by  a deconvolution  of  each 
of  the  peaks  of  the  PFD.  If  we  can  deconvolute  each  of  the 
peaks,  the  distribution  of  interatomic  distances  in  vitreous 
silica  can  be  directly  obtained.  The  distribution  of  atomic 
distances  in  glass  is  described  by  the  Distance  Distribution 
Functions  (DDF) . 
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Pete rmination  of  Vitreous  S ilica  DDF 

Theoretically,  the  DDF  can  be  determined  by  deconvolv- 
ing each  of  the  peaks  of  the  vitreous  silica  PFD.  However, 
this  procedure  is  not  practical  because  of  the  overlapping 
of  the  PFD  peaks.  In  practice,  the  following  equivalent 
procedure  was  adopted  for  determining  the  vitreous  silica 
DDF.  First,  a set  of  relevant  pair  functions  [Warren 
C 1 9 6 9 a , b ) ] , i.e.,  those  for  Si-0,  0-0,  and  Si-Si  interac- 
tions were  calculated.  Assuming  various  DDF,  the  DDF  were 
convoluted  with  the  corresponding  pair  functions.  The  con- 
voluted pair  functions  together  with  the  proper  lb  ^ and  r^^ 
values  were  used  to  calculate  the  left-hand  side  of  Eq.  (3.2). 
The  set  of  DDF  used  was  varied  until  a good  fit  was  obtained 
between  the  experimental  and  the  calculated  PFD.  A good 
fit  could  be  obtained  only  up  to  the  fourth  peak  of  the  PFD. 

A typical  vitreous  silica  DDF  obtained  by  Mozzi  and 
Warren  is  shown  in  Fig.  3.  Figure  3a  gives  the  DDF  of  Si- 
1 Si  atoms,  while  the  corresponding  distributions  of  the 
Si-O-Si  bond  angles,  V(a),  is  shown  in  Fig.  3b.  The  ratio 
of  the  0-lst  0 to  the  Si-lst  0 distance  is  equal  to  /WJ 3 
and  is  exactly  what  would  be  expected  for  tetrahedral  bond- 
ing. Hence  the  unit  of  structure  in  vitreous  silica  can 
be  taken  as  a SiO^  tetrahedron,  with  silicon  at  the  center 
and  four  oxygens  at  the  corners  of  the  tetrahedron.  The 
Si-ls1  Si  distance  is  the  distance  between  the  Si  atoms  at 
the  centers  of  two  adjacent  tetrahedra  (Fig.  4).  It  can  be 
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Figure  3.  Functions  describing  Si-1  Si  distance  distribution  in 
vitreous  silica  [adapted  from  Mozzi  and  Warren  (1969)]. 
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Figure  4.  Schematic  of  two  adjacent  tetrahedra  in 
vitreous  silica  showing  the  Si-lst  Si 
distance . 
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s "t 

shown  that  the  Si-1'  Si  DDF  is  related  to  the  distribution 
V(n)  of  Si -0-Si  bond  angle  a by 

V (a)  « G(rii)  1.62  cos  a/2  (3.7) 

The  DDF  of  the  third  and  subsequent  peaks  are  dependent  on 
one  another  because  of  the  dependence  of  the  corresponding 
bond  lengths  on  a. 

Once  a set  of  DDF  giving  a calculated  PFD  in  close 
agreement  with  the  experimental  PFD  is  obtained,  we  have  a 
quantitative  statistical  description  of  the  structure  of 
vitreous  silica.  The  DDF  are  essentially  similar  to  the 
atomic  density  distribution  functions  obtained  by  the  old 
RDF  methods.  But  because  the  PFD  analysis  is  more  rigorous, 
the  DDF  give  a more  precise  description  of  the  structure  of 
vitreous  silica.  We  could  describe  a particular  structural 
state  of  silica  by  a set  of  DDF  or  by  a characteristic 
feature  of  the  DDF  such  as  the  half-maximum  breadth  of  DDF 
curves.  Changes  in  the  structure  of  vitreous  silica  can  be 
described  by  the  corresponding  changes  in  the  DDF.  This  is 
the  approach  used  in  the  present  investigation  to  study  the 
structural  changes  in  vitreous  silica  as  a function  of  heat- 
treatment,  water  content  and  neutron  irradiation. 


CHAPTER  IV 


EXPERIMENTAL  PROCEDURE 
Materials 

Various  grades  of  commercial  silica  were  used  during 
the  course  of  the  present  investigation.  The  details  in- 
cluding the  method  of  manufacture,  impurity  and  water  con- 
tent are  given  in  Table  2. 


Heat -Treatment  Procedures 

The  heat-treatments  were  carried  out  on  samples  of 
"T-D"  clear  fused  quartz.  The  initial  heat-treatments 
were  confined  to  bulk  samples.  No  change  was  detectable 
in  the  x-ray  diffraction  patterns  of  the  bulk  samples  even, 
after  long  periods  of  heat -treatment  at  various  tempera- 
tures between  1400°  and  166 1°C.  Hence  it  was  decided  to 
try  heat  - treatment s of  powder  samples  of  "T-D"  clear  fused 
quartz.  The  powder  diffraction  patterns  began  to  change 
after  a definite  period  of  time.  The  details  of  the  heat- 
treatment  procedures  are  discussed  in  this  section. 
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Details  of  the  Materials  Used  in  the  Present  Investigation 


28 


+-> 

o\o 

P 3 

to 

CD  CD 

o\© 

o 

o\o 

o\« 

P P 

OvJ 

o 

H* 

rsj 

3 3 

rH 

o 

o 

rH 

©NP 

E:  o 

• 

• 

• 

<N2 

a 

O 

o 

o 

O 

rH 

LO  LO 

5 

3 

4 

o H-  3 

O 

3 

E 

to  3 

rH 

3 

E 

O PL, 

Pi 

.3 

3 PD 

•H  CD 

P 

P 

CD  PD 

rH 

•H 

r-H 

K)  K) 

to 

to 

+->  rH 

3 to 

o o 

CM  O 

to 

to 

•H 

P 0) 

CM  CMO  O CM 

CD 

3 

^ r: 

a>  -h 

rH  (L) 

•H  NtS 

rH  3 rH 

rH 

rH 

3 ctf 

p 

< ch 

h^2; 

< 2 

< 

•H 

to 

to 

S 4-J 

P P 

t/)  </) 

to  to  to 

P 

s 

p 

E 

•H 

O 3 

Cl  Gj 

3 3 3 

CD 

Cl, 

CD 

PD 

</) 

•!—,  Ph 

PS 

PL, 

PS 

PD 

rH  (/) 

3 e 

rH  0 

• H 3 

P 

P 

rH  0 

s *-• 

< PD 

T 

K 

N 

O rH 

O rH 

rH 

CD 

P 

1 

1 

PS 

I 

H 

to 

to 

o 

p 

X 

X 

o 

rH  t 

u 

p 

p 

p 

O 3 

3 

o 

CD  O 

p 

<P 

O 

tp 

33  3 

3 

PS 

3 p* 

o 

Sh-H 

3 

u 

3 O 

»*3 

3 

o 

P O 

3 

H-  t 

S 

p 

3 P 

o 

CD  rH 

Mh 

•H 

to  o 

'-H 

<p 

<p 

to 

3 -H 

O 

o 

o o 

3 

PS  r/j 

•H 

m rH 

PD 

33 

3 

P 3 

3 

m • 

o 

o 

P O 

CD 

P 

P O CD 

PS 

•H 

U -H 

E 

to 

o E 

P 

to  rH 

!D  3 rH 

3 

Sh 

CD  to  3 

CD 

3 3 

rH  3 3 

rH 

P 

■3  »H  i — 1 

S 

PD  P 

w m P 

IP 

U 

>■  to  <P 

to 

CD 

•H 

P 

•H 

J-i 

3 

P- 

e 

•H 


3 

+-> 

O 

H 


33 

3 

P> 

Pi 

CD 

N ■ 

rH 

PQ 

rd 

to 

P P 

•H 

i 

•H 

PD 

3 P 

to 

Pi 

, — l 

o 

CD  3 

O 

r — t 

•H 

Gj 

CD 

rH  3 

CD 

•H 

to 

i — 1 

o cr 

P 

to 

O 

u 

PD 

P 

o 

P 

•H 

E 

P 33 

•H 

CD 

P 

U 

3 

Q <D 

> 

P 

o 

•H 

CO 

1 to 

P 

CD 

rH 

E-i  3 

Pi 

•H 

CD 

•H 

r <p 

KH 

> 

CO 

CO 

>> 

P to 

P 

33 

CD  P 

CD 

3 

h3  3 

P 

P 

£ CD 

3 

CO 

O E 

S 

PDP 

<P 

CO 

<P 

o 

HD  CD 

O 

3 p 

P 

tD 

3 P 

P 3 

P 

i 

O CD 

3 

PS  P 

CD  4-* 

p 

rH  3 

DH  3 

3 

3 CD 

<P  O 

SP 

M PS 

W O 

less  than  0.04 


Table  2 (con 


29 


t — \ 

p 


p 

p p 

0 0 

1 

p p 

1 

p p 

1 

^ o 

u 

E 

U Ph 

•H  pi 

rH 

rH  •> 

P P 

P 0 

1 

0 -H 

1 

* 2 P 

1 

•H 

P P 

O P 

■n  pi. 

p e 

2 w 

0 

p 

p 

p 

u 

p 

4-1 

P 

P 

P 

2 

4H 

1 

o 

1 

M3 

o 

rP 

P 

0 

2 

m3 

0 

P 

P> 

0 

o 

i — 1 

rH 

Pi 

•H 

£ 

V) 

p 

O 

CO 

•P 

P< 

O 

i 

m3 

•H 

P 

P 

nj 

CO 

4-i 

hH 

o 

£ 

0 

O 

P 

^ £ 

P 

P o 

P 

3 -H 

P 

0 P 

2 

03 

< 


co 


2 

o 

p 

p 

0) 

•n> 

£ 

(D 


m3 

rH 

0 


P 

•H 

Pi 


O 

Pi 

t— i 


N 

P 

P 

ci 

P 

cy 

p 

CD 

P 

p 

Pi 

p 

0 

P 

P 

p 

co 

m3 

0 

P 

•H 

P 

O 

P 


X 

P X 
•H  X> 

£ 

P 

~ P £ 

NO® 

6Hp 

u 

\ X - 

O f — I i — I 
N P H 
O p 
i — I 

<+■1 


■H  K 


2 
O P P 


P 

P 


0 

^ I — I 

P Pi  P 
H E ^ 
4-1  P 
CO  *> 


i — i co 

P 0 

P • -H 
O > P 
P 0 o 
« p 

p p 

O P 

O rH  O 

P -O 

• P P 

M3  p P PI 

P O 43 

P -H  P 0 

<H  P P 

bO  P P O 

P 'HOP 

W 3P  ft 
P P 0 
i P 0 rH 

T3  P P 0 

p -H  bOH 

0 

P P P rH 

I — i 0 0 t — t 

rH  P H 0 

p <p  tiom 

^ p p 

0 4-1 


M3 

P 

o 

0 

P 

0 

P 

P 

0 

P 

b® 

t/) 

U 

0 

P 

3 

• 

•H 

P 

•H 

oj 

< 

M3 

O 

> 

• 

P 

P 

p 

CO 

0 

rP 

• 

CO 

P 

• 

p> 

P 

H 

1 — 1 

P 

P 

p 

0 

o 

• 

E 

i — i 

p 

•""3 

0 

P 

P, 

p 

p 

0 

E 

p 

• 

p 

rP 

p 

0 

P 

0 

H 

CO 

p 

Q 

O 

P o 


30 


Bulk  Heat -Treatments 

The  bulk  samples  were  heat-treated  at  a series  of  tem- 
peratures between  1400°  and  1661°C  for  various  times  ranging 
from  29  to  360  hours.  The  sample  thickness  varied  from 
1-1.5  cm.  The  heat -treatment  details  are  summarized  in 
Table  3.  In  order  to  reduce  surface  crystallization  at 
temperatures  above  1500 °C,  the  samples  were  flame-polished 
and  heated  in  an  amorphous  carbon  crucible  in  an  argon  at- 
mosphere graphite  tube  furnace.  The  samples  were  inserted 
into  the  quench  chamber,  the  furnace  flushed  with  aigon, 
and  the  samples  inserted  into  the  hot  zone  at  temperature. 
After  the  desired  heating  interval,  the  samples  were  with- 
drawn into  the  quench  chamber  and  cooled  overnight.  Sam- 
ples heat-treated  below  1500 °C  were  placed  in  Pt  boats  and 
heated  in  a SiC  muffle  furnace  in  ambient  atmospheres. 

Sample  surface  preparation  for  x-ray  diffraction 

The  heat -treatment s invariably  produced  surface  crys- 
tallization. The  crystallized  surface  was  removed  and  the 
surface  polished  before  running  the  diffraction  pattern. 

The  effect  of  polishing  on  the  diffraction  data  was  deter- 
mined as  follows.  A bulk  sample  heat-treated  for  40  hours 
at  1450°C  was  taken  through  the  following  polishing  se- 
quence: 120  grit  SiC  powder  on  a polishing  wheel,  320, 

400  and  600  grit  SiC  paper,  0.3y  alumina  and,  finally,  an 
optical  polish  with  ceria.  The  first  peak  in  the  Cu 
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Table  3 

Heat -Treatment  Procedures  for  Bulk  Samples 


Temperature 

(°C) 

Time 

(hours) 

Furnace  Used 

Atmosphere 

1 4 0 0 a 

40-120 

SiC  muffle^ 

Ambient 

1450a 

40-360 

SiC  muffle 

Ambient 

1595° 

58-293 

Graphite  resis- 
tance tube^ 

Commercial  purity 
argon 

4 

1 6 6 1 c 

29 

Graphite  resis- 
tance tube 

Commercial  purity 
argon 

Samples 

placed  in 

Pt  boats. 

l. 

DPereny  Corporation. 

Samples  flame -polished  and  supported  in  an 
amorphous  carbon  crucible. 

^Astro  Industries,  Inc. 
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diffraction  pattern  was  scanned  after  each  of  the  above 
polishing  steps.  It  was  found  that  the  first  peak  was 
essentially  unchanged.  The  first  peak  diffraction  data 
after  the  120  grit  SiC  polish  were  compared  with  the  data 
obtained  after  the  optical  polish  (Fig.  5).  It  was  seen 
that  the  intensity  was  essentially  unaltered  as  a result  of 
polishing.  The  data  for  the  intermediate  steps  were  found 
to  be  essentia.lly  the  same.  The  sample  was  also  given  an 
etching  treatment  with  an  HF  solution  (50  percent  IIF, 
remainder  HNO_)  for  one  minute.  The  etching  step  had  no 
effect  on  the  first  peak  intensity  data.  The  polishing  se- 
quence was  stopped  at  the  0.3p  alumina  step  for  the  other 
bulk  heat-treated  samples. 

Powder  Heat-Treatments 

"T-D"  clear  fused  quartz  samples  were  crushed  with  a 
hammer,  boiled  overnight  in  dilute  HC1  to  dissolve  the  iron, 
washed,  dried,  and  ground  in  an  agate  mortar  to  pass  through 
a 325-mesh  screen.  The  fine  (-325  mesh)  powders  were  washed 
with  triply  distilled  water  to  minimize  surface  contami- 
nants, and  dried  overnight  in  an  oven  at  200°C.  The  x-ray 
diffraction  pattern  of  the  powder  sample  was  identical  to 
that  of  the  bulk  sample  of  "T-D"  clear  fused  quartz.  The 
powders  x^ere  mounted  in  a Pt  boat  and  heat-treated  in  a SiC 
muffle  furnace  at  1200°C  (±5°C).  The  heat -treatments  were 
carried  out  in  ambient  atmospheres.  No  change  was 
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Figure  5.  Effect  of  polishing  on  the  first  diffraction 
peak  of  vitreous  silica  heat-treated  40  hours 
at  14  50  °C . 
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observable  in  the  diffraction  patterns  after  14  hours  at 
temperature;  a slight  change  was  a noticeable  after  about 
19  hours  and  the  change  became  pronounced  after  30-1/2  hours 
at  temperature.  A PFD  analysis  was  performed  on  the  powders 
heat-treated  for  19  and  30-1/2  hours. 

The  powder  samples  for  x-ray  diffraction  were  made  by 
packing  the  powders  in  a specially  designed  sample  holder 
made  of  Teflon.*  The  sample  holder  was  designed  with  a 
rectangular  cavity  such  that  when  the  silica  powders  were 
packed  in  the  cavity,  there  was  no  x-ray  scattering  from 
either  the  sides  or  the  back  of  the  cavity.  The  particle 
size  distribution  in  the  powders  was  such  that  the  powders 
could  be  packed  in  the  holder  without  any  binder.  The 
powders  tended  to  agglomerate  on  heat -treatment , but  this 
could  be  corrected  by  a slight  grinding,  before  packing 
into  the  holder  for  x-ray  diffraction  studies. 


X-ray  Diffraction  Procedures 

In  their  investigation  of  the  structure  of  vitreous 
silica,  Mozzi  and  Warren  used  primary  beam  monochromated 
Cu  and  Rh  radiations.  They  used  the  Warren-Mavel  (1965) 
fluorescence  technique  to  experimentally  eliminate  the 
Compton  scattering  when  using  Rh  radiation.  During  the 


^Trademark  of  E.  I.  DuPont  de  Nemours,  Inc. 


initial  stages  of  the  present  investigation  an  attempt  was 
made  to  duplicate  the  system  used  by  Mozzi  and  Warren. 

The  fluorescence  technique  did  not  produce  sufficient 
intensity  due  to  the  lack  of  a suitable  LiF  crystal  mono- 
chromator. Therefore  a standard  Norelco  goniometer  with  a 
diffracted  beam  monochromator  and  Mo  radiation  was  modified 
to  obtain  the  experimental  intensity  values  for  high  S 
values.  The  modification  helped  to  eliminate  Compton  scat- 
tering at  high  angles. 

Copper  Diffraction  System 

The  diffraction  patterns  at  low  S values  (0.11-1.2  A ^ 
were  obtained  using  primary  beam  monochromated  Cu  radiation 
A GE  x-ray  tube  and  spectrogoniometer , together  with  a 
Norelco  Xe  proportional  counter  and  Norelco  electronics 
system  consisting  of  a linear  amplifier,  pulse  height  analy 
zer  and  scaler-ratemeter , were  used  for  obtaining  the  coppe 
diffraction  data.  The  diffraction  patterns  were  scanned 
from  150°  to  10°  20  at  a rate  of  2°  per  minute.  The  dif- 
fracted intensity  was  recorded  using  a Bristol  chart  re- 
corder with  a chart  speed  of  1/2"  per  minute.  Four  scans 
were  made  for  each  sample  and  the  average  of  the  four  scans 
determined.  The  average  intensity  was  then  corrected  f Gl- 
air-scattering , beam-off  background,  polarization  factor 
and  normalized  using  a modification  [Warren  (1969a,b)]  of 
a normalization  procedure  originally  proposed  by  Norman 
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(1957).  The  air-scattering  and  background  corrections 
were  determined  as  follows.  A scan  was  made  with  the  sam- 
ple holder  removed  but  the  x-rays  on.  The  intensity  re- 
corded is  due  to  the  background  noise  in  the  system  and 
scattering  from  a volume  of  air  twice  that  present  when  the 
sample  is  in  position.  The  background  noise  can  be  evalu- 
ated by  shutting  off  the  x-rays.  Corrections  for  the  in- 
strumental broadening  effect  can  be  carried  out  by  a modi- 
fied Sauder  method  [Warren  and  Mozzi  (1970)].  The  Sauder 
method  applied  to  the  Cu  data  from  vitreous  silica  revealed 
that  the  instrumental  broadening  was  negligible.  Therefore 
the  average  intensity  had  to  be  corrected  only  for  air- 
scattering  and  background  noise.  The  corrected  average 
intensity  was  then  divided  by  the  polarization  factor,  p.f., 
given  by 

1 + cos^  20.,  cos^  20 

p.f.  = H (4.1) 

1 + cos^  20 

where  20^  is  the  monochromator  angle.  The  corrected  inten- 
sity was  normalized  by  the  Norman  method.  The  details  of 
the  normalization  procedure  are  presented  in  Appendix  1. 

Molybdenum  Diffraction  System 

The  diffraction  data  for  values  of  S greater  than 
1.2  X"1  were  obtained  using  MoKa  radiation.  The  experi- 
mental system  consists  of  a standard  Norelco  vertical 


37 


goniometer  with  a LiF  curved  crystal  diffracted  beam  mono- 
chromator. A Norelco  scintillation  counter  was  used  to 
monitor  the  diffracted  intensity.  The  pulse  height  analyzer 
was  set  to  allow  only  90  percent  transmission  of  the  MoKa 
radiation . 

The  diffracted  intensity  was  monitored  using  a fixed 
count  procedure.  The  time  needed  for  1,000  counts  was 
measured  as  a function  of  angle,  in  2°  steps  from  150°  to 
80°  20  and  in  1°  steps  from  80°  to  42°  26.  The  relative 
probable  counting  error  was  thus  independent  of  angle.  The 
measured  intensity  was  corrected  for  air-scattering,  beam- 
off  background,  polarization  factor  and  multiple  scattering 
and  then  normalized.  The  Compton  scattering  was  negligible 
in  the  S value  range  of  interest.  The  air-scattering  and 
beam-off  background  contributions  were  determined  in  a man- 
ner similar  to  that  for  the  Cu  diffraction  data.  The  polar- 
ization factor  (p.f.)  for  the  diffracted  beam  monochromator 
system  is  given  by 

1 + cos^  20w  cos^  20 

p.f.  = S (4.2) 

1 + c.os^  2 6^ 

Multiple  scattering  is  significant  for  the  case  of  Mo  radi- 
ation scattered  by  vitreous  silica.  The  multiple  scatter- 
ing correction  was  evaluated  using  the  generalized  equation 
derived  by  Warren  and  Mozzi  (1966)  for  the  ratio 
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1(2)  _ B2  Q C 2 6 ,q  ,b) 
1(1)  " J(P)  E ATyT (m) 


(4.3) 


where 

1(2)  and  I 

(1) 

refer 

to  the  doubly 

and  singly 

seat- 

tercd 

intensity , 

B = 

E 

Z^2,  and  Q is 

a 

function 

ed  f 20 

u . c . 

q and 

b.  Warren 

and 

Mozzi 

have  developed 

a relation  be- 

tween 

J(p) , qs  0 

and 

"b". 

J(p)  = B[q  + (1  -q)  / (1-b  sin26)]  (4.4) 


where  ,fp!i  is  a function  of  the  scattering  angle  20.  J(p)  can 
be  assumed  to  be  the  same  as  the  independent  unmodified  scat- 
tering Z f.z'  as  the  Compton  is  effectively  eliminated, 
u . c . - 

"q"  and  "b"  can  be  calculated  from  the  Eq.  (4.4).  Knowing 
"q"  and  "b",  "Q"  can  be  determined  from  the  table  given  by 
Warren  and  Mozzi  (1966)  , (Table  4) . and  y^ (m)  are  the 

atomic  weight  and  mass  absorption  coefficient  of  atom  type 
"i".  The  values  of  the  ratio  I(2)/I(l)  were  calculated  for 
20  values  of  60°,  90°,  120°  and  150°.  These  values  (Fig.  6) 
were  plotted  against  20,  and  a smooth  curve  drawn  connecting 
the  calculated  points.  The  values  of  the  ratio  of  the  doubly 
scattered  to  the  singly  scattered  beam  were  read  off  from 
the  graph  as  a function  of  20  in  2°  steps  from  150°  to  80° 
and  1°  steps  from  80°  down.  The  correction  for  multiple 
scattering  could  then  be  easily  achieved  by  dividing  the 
corrected  intensity  values  by  the  corresponding  values  of 
the  ratios  I (1)  + I (2)/I (2) . It  should  be  noted  that  only 
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Figure  6.  Ratio  of  double  to  single  scattering  contribu- 
tions from  vitreous  silica  (MoKa) . 


29 

30 

60 

90 

120 

150 

180 

30 

60 

90 

120 

150 

180 


40 


Table  4 

Values  of  Q(20 ,q ,b)xl0^a 


b = 10  b = 20  b = 40  b = 60  b = 80  b = 100 


q = 0.00 


192 

104 

50.8 

196 

87 

33.7 

164 

64 

22.2 

124 

45 

14.9 

10  3 

36 

11.9 

98 

34 

11.1 

q = 0.05 


219 

126 

68.1 

238 

121 

59.1 

218 

105 

51.9 

179 

86 

43.9 

156 

75 

39.4 

150 

73 

38.1 

31.3 

21.6 

16.0 

18 . 3 

11.7 

8.1 

11.5 

7.1 

4.8 

7.6 

4.6 

3.2 

6.0 

3.7 

2.5 

5.6 

3.4 

2.3 

46.2 

34.9 

2 8.1 

39.4 

30.1 

24.7 

35.7 

28.1 

23.7 

31.1 

25.0 

21.4 

28.3 

23.0 

19.8 

27.5 

22.4 

19.4 

aAdapt.ed  from  Warren  and  Mozzi  (1966). 
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the  doubly  scattered  component  is  taken  into  account.  The 
higher  order  scattering  is  generally  negligible.  The  Mo 
data  were  then  normalized  as  follows. 

Strong  and  Kaplow  (1966)  have  shown  that  for  I^O, 
glass  the  Compton  scattering  can  be  essentially  eliminated 
by  using  a diffracted  beam  monochromator  with  MoK  radia- 
tion. In  the  present  investigation  a preliminary  study  was 
made  to  determine  the  effectiveness  of  various  slit  systems 
in  eliminating  Compton  scattering.  Using  a 1°  divergence 
slit  and  a 4°  receiving  slit,  the  Mo  diffraction  patterns 
were  recorded  between  90°  and  150°  28,  using  1°,  1/2°, 
0.006"  and  0.003"  scatter  slits.  The  results  are  summar- 
ized in  Fig.  7.  It  can  be  seen  that  a 0.006"  scatter  slit 
effectively  brings  down  the  tail  portion  of  the  Mo  diffrac- 
tion pattern,  without  much  loss  in  intensity.  As  Compton 
scattering  is  mainly  responsible  for  the  high  angle  por- 
tion of  a diffraction  pattern  going  up,  it  can  be  seen  that 
the  l°4°0.OO6"  slit  system  effectively  eliminated  the  Comp- 
ton scattering.  The  actual  experiments  were  performed 
using  a Mo  x-ray  tube  run  at  45  KV,  18  mA  setting,  with 
the  1°4°0.006"  slit  combination.  This  resulted  in  satis- 
factory elimination  of  Compton  scattering. 

The  efficiency  of  removal  of  Compton  scattering  by 
this  method  was  then  determined  as  follows.  A constant 
was  determined  such  that  the  corrected  molybdenum  data 
multiplied  by  the  constant  coincided  with  the  normalized 
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Figure  7. 


Effect  of  scatter-slit  size  on  the  high  angle 
region  of  the  Mo  diffraction  pattern. 
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total  intensity  curve  determined  using  Cu  radiation  at  S 
values  near  zero.  However,  the  molybdenum  data  curve  be- 
gins to  diverge  from  the  copper  curve  at  higher  values  of 
S due  to  the  progressive  elimination  of  the  Compton  scatter- 
ing at  these  S values.  The  difference  between  the  normal- 
ized copper  data  and  the  molybdenum  data  was  assumed  to 
be  entirely  due  to  the  Compton  scattering  which  has  been 
eliminated.  This  difference  was  measured  as  a function  of 
S.  Now,  the  total  Compton  scattering  from  vitreous  silica 
can  be  obtained  by  the  use  of  the  International  Tables  for 
X-ray  Crystallography,  Vol.  Ill  (1962)  and  hence  the  frac- 
tion of  Compton  scattering  which  is  removed  can  be  deter- 
mined as  a function  of  S.  The  fraction  of  Compton  scatter- 
ing which  remains  is  of  interest  to  us  and  is  shown  in 
Fig.  8.  The  data  were  not  sufficient  to  cover  the  whole 
region  of  molybdenum  S values  and  so  a rather  extensive 
extrapolation  was  necessary.  The  normalized  molybdenum 
data  were  then  corrected  for  the  fraction  of  Compton  scat- 
tering remaining,  again  using  the  total  values  calculated 
from  the  tables. 

Apart  from  the  assumption  of  complete  coincidence  of 
normalized  copper  and  molybdenum  diffraction  data  at  low 
S values,  this  procedure  suffers  from  the  following  draw- 
backs as  well. 

1.  The  extrapolation  of  the  Compton  scattering  cor- 
rection data  is  made  essentially  over  the  entire  molybdenum 
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Figure  8. 


Ratio  of  Compton  remaining  (ip)  to  total 
Compton  (iT)  for  1°4°0.006"  slit  system. 
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diffraction  region  and  the  accuracy  of  this  extrapolation 


is  suspect. 


2.  A great  reliance  is  placed  on  the  tabulated 
Compton  scattering  values.  Alexandropoulos  (1966)  experi- 
mentally determined  the  Compton  scattering  from  a vitreous 
silica  sample  at  an  S value  of  2 using  molybdenum  radia- 
tion. He  found  that  the  experimental  Compton  scattering 
value  was  less  than  the  value  calculated  from  the  tables 
by  about  15  percent.  Thus  an  unknown  error  is  introduced 
in  the  extrapolation  procedure. 

3.  The  differences  between  the  normalized  Cu  and  Mo 
curves  were  fluctuating  at  S values  where  the  diffraction 
curve  undergoes  rapid  changes  resulting  in  some  scatter  in 
the  data  points  of  Fig.  8. 

Because  of  the  above-mentioned  difficulties,  the  follow- 
ing alternative  procedure  was  tried.  The  extrapolation  pro- 
cedure indicates  that  the  fraction  of  Compton  scattering 
which  remains  (iR/i'T)  in  the  region  of  S values  above  1.2  A 1 
is  a very  small  fraction  of  the  total  Compton  in  this  re- 
gion. The  assumption  could  therefore  be  made  that  the  Comp- 
ton scattering  was  completely  eliminated  in  this  region  and 

the  corrected  Mo  intensity  data  were  normalized  to  oscillate 

2 

about  the  independent  unmodified  scattering  curve  £ f. 


u . c . 


1 


in  this  regaon.  This  procedure  gave  rise  to  a normalized 
Mo  data  curve  which  gave  a good  fit  to  the  normalized  un- 
modified data  for  Cu  in  the  region  around  S values  of 
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1.2  A 1 (Fig.  9).  The  error  in  neglecting  the  fraction  of 
Compton  remaining  at  values  of  S greater  than  1.2  A 1 could 
be  taken  care  of  by  the  correction  procedure  outlined  in 
the  next  section.  The  alternative  procedure  was  used  for 
all  the  samples  studied  in  the  present  investigation  except 
for  "T-D"  clear  fused  quartz,  for  which  the  original  pro- 
cedure was  used. 


PFD  Calculations 

The  normalized  intensity  data  in  electron  units  per 
unit  of  composition  as  a function  of  k(2iTS)  were  used  to  cal- 
culate the  integral  in  the  PFD  equation 

N.  . ? 

X X — ^ P • • (r)  = 2tt  rp  X Z. 
u.c.  i.  lj  J u.c.  -1 

k 2 2 

+ /Qm  ki(k)e  a ^ sin  rk  dk  (3.2) 

The  Fourier  integral  of  Eq . (3.2)  was  evaluated  using  a 
modification  of  a computer  program  developed  by  Dr.  Heritage 
(1969)  at  the  University  of  Florida.  The  average  electron 
density  in  the  sample  depends  on  the  overall  density  of  the 
sample  under  consideration.  The  average  electron  density 
was  taken  to  be  the  same  for  all  the  samples  except  for 
irradiated  silica  and  silicic  acid  which  have  densities 
significantly  different  from  the  other  samples.  The  com- 
puter program  was  written  to  produce  a printout  and  plot  of 
PFD  as  a function  of  r in  X. 
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Figure  9.  Normalized  copper  (unmodified)  and  Mo  intensity  data  showing  mat 
at  S values  around  1.2  A"l. 
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The  computer  output  plot  which  is  a plot  of  the  experi- 
mental PFD  in  general  contains  ripples  especially  in  the 
region  of  "r"  values  below  the  first  peak.  These  ripples 
are  caused  by: 

1.  errors  in  the  experimental  intensity  values, 

2.  errors  in  the  normalization  of  the  experimental 
intensity  values, 

3.  errors  in  the  tabulated  atomic  scattering  factors 
and  Compton  scattering  factors,  and 

4.  termination  of  the  Fourier  transform  at  a finite 
"k"  value. 

The  ripples  can  cause  some  difficulty  in  the  interpretation 
of  the  experimental  PFD.  The  ripples  due  to  the  errors  men- 
tioned above  can  be  eliminated  by  a procedure  originally 
proposed  by  Kaplow,  Strong  and  Averbach  (1965).  Appendix  2 
gives  an  outline  of  the  correction  procedure.  The  termina- 
tion ripples  are  still  present  after  application  of  this 
procedure  but  they  cause  little  harm  because  the  same  rip- 
ples are  present  in  the  calculated  PFD  also.  Experimental 
PFD  of  vitreous  silica  before  and  after  the  correction  pro- 
cedure (to  be  called  Kaplow  correction  procedure)  are  shown 
in  Fig.  10.  It  is  clear  that  the  correction  procedure  has 
eliminated  most  of  the  ripples  without  significantly  alter- 
ing the  position  or  shape  of  the  main  experimental  peaks. 

The  Kaplowr  correction  procedure  was  used  to  eliminate  the 
error  ripples  for  all  the  PFD  determined  in  the  present 
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investigation.  A comparison  of  the  PFD  before  and  after 
the  correction  procedure  was  made  in  each  case  and  con- 
firmed that  the  correction  procedure  does  not  alter  the 
position  or  shape  of  the  PFD  peaks.  A slight  change  in  peak 
values  due  to  redistribution  of  the  ripples  was  the  only 
significant  change. 


CHAPTER  V 


RESULTS  AND  DISCUSSION 

The  results  of  the  present  investigation  are  presented 
and  discussed  in  this  chapter.  The  effect  of  heat-treat- 
ment on  the  x-ray  diffraction  patterns  of  the  bulk  samples 
is  discussed  first,  followed  by  an  analysis  of  the  diffrac- 
tion data  of  the  heat-treated  powder  samples.  The  PFD  of 
the  samples  with  different  water  contents  are  then  presented 
and  discussed.  A DDF  analysis  of  the  structural  changes  in- 
duced in  a sample  of  Optosil  due  to  heavy  neutron  irradia- 
tion completes  this  chapter. 


Heat -Treatment  Results 


The  diffraction  patterns  of  the  bulk  samples  were  not 
changed  by  the  heat -treatments  used  in  the  present  investi- 
gation. A comparison  of  the  normalized  intensity  data  for 
a bulk  sample  of  "T-D"  clear  fused  quartz  before  and  after 
a heat-treatment  of  360  hours  at  1450°C  is  given  in  Fig.  11. 
It  is  clear  that  there  is  no  significant  difference  between 
the  two  curves.  Internal  crystal  growth  was  observed  in 
some  of  the  heat-treated  samples  (Table  5)  but  no 
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Figure  11.  Normalized  intensity  curves  for  "T-D"  clear  fused  quartz  and 
bulk  sample  heat-treated  for  360  hours  at  1450°C. 
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Table  5 

Growth  of  Internal  Crystals  in 
Heat-Treated  Bulk  Samples 


Heat-Treatment  Internal  Crystals 

40-120  hours  at  1400°C  No  internal  crystals  observable 

40-360  hours  at  1450°C  No  internal  crystals  observable 

58-293  hours  at  1595°C  Internal  growth  of  crystals 

observed 


29  hours  at  1661°C 


Internal  crystals  present 
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significant  changes  in  their  diffraction  patterns  could  be 
detected , 

X-ray  diffraction  patterns  of  the  powdered  vitreous 
silica  samples,  on  the  other  hand,  show  considerable 
changes  as  a function  of  time  at  the  heat -treatment  temper- 
ature. The  diffraction  patterns  begin  to  show  a change 
after  about  19  hours  at  1200°C  (±5°C),  and  the  changes  be- 
come more  pronounced  after  30-1/2  hours  at  temperature. 

The  normalized  intensity  curves  for  "T-D"  clear  fused  quartz 
(0  hours)  and  powder  samples  heat-treated  for  19  and  30-1/2 

hours  at  1200°C  (±5°C)  are  shown  in  Fig.  12.  With  increas- 

° - 1 

ing  time  of  heat -treatment , the  first  peak  at  S ~ 0.25  A 
becomes  narrower  and  increases  in  intensity.  Clearly  de- 
finable Bragg  peaks  at  positions  corresponding  to  low  cris- 
tobalite  begin  to  appear  after  30-1/2  hours. 

The  Kaplow  corrected  PFD  of  the  three  samples  computed 
using  the  right-hand  side  of  Eq . (3.2)  are  shown  in  Fig.  13. 
Significant  changes'  are  noticed  in  the  peaks  of  the  PFD  with 
increasing  time  of  heat -treatment . Due  to  the  nature  of  the 
Kaplow  correction  procedure  (Appendix  2) , the  cut-off  value 
S of  the  samples  PFD  varies  and  a direct  comparison  of 
the  PFD  is  difficult.  However,  the  DDF  approach  takes  the 

variation  of  S values  into  account  and  hence  is  a more 
max 

reliable  index  of  structural  changes.  Hence  the  PFD  were 
analyzed  by  the  DDF  approach.  The  PFD  of  "T-D"  clear  fused 
quartz  has  the  same  peak  positions  as  the  PFD  obtained  by 
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Figure  12.  Normalized  intensity  curves  for  "T-D"  clear  fused  quartz  (0  hours) 
and  powder  samples  heat-treated  for  19  hours  and  30-1/2  hours 
at  1200  °C . 


Figure  13.  Experimental  PFB  (Kaplow  corrected)  of  "T-B"  clear  fused 
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Mozzi  and  Warren  for  their  sample  of  vitreous  silica.  How- 
ever, the  magnitudes  of  the  peak  values  vary.  Hence  a com- 
parison could  be  made  only  after  obtaining  the  DDF  of  the 
"T-D"  cleai*  fused  quartz  PFD. 

Determination  of  the  DDF  for  Heat-Treated  Samples 

The  first  peak  of  the  PFD  of  the  three  samples  all  lie 
at  about  1.62  A.  Assuming  tetrahedral  coordination  of  the 
oxygen  atoms  about  the  silicon  atoms,  the  position  of  the 
second  and  subsequent  peaks  can  be  calculated.  The  DDF  of 
the  first  two  peaks  were  found  to  be  the  same  as  used  by 
Mozzi  and  Warren.  The  DDF  of  the  third  and  subsequent 
peaks  were  calculated  as  follows.  The  Si-lst  Si,  Si-2nd  0, 
Si-2nd  Si  and  0-2nd  0 distances  in  vitreous  silica  can  be 
described  by  the  following  equations: 


Si - 1st  Si 

r . 
id 

= 3.24  sin 

a/ 2 

(5.1) 

Si - 2nd  0 

? 

rt . 
13 

= 9.65-7.00 

cos 

a - 4.96  sin  a cos  B (5.2) 

0- 2nd  0 

2 

r-  . 

= 4.654(cos 

a - 

sin  B sin  y + cos  a cos  B cos  y) 

ij 

+ 7.015  [2-2  cos  a - 0.9417  sin  a(cos  B + cos  y)  ] 

(5.3) 

Si  -2nd  Si  t2.,  = 10 . 5 (sin  a/2 ) 2 [ 2 . 673- 0 . 94 17  sin  a(cos  B 
J 

+ cos  y)]  - 5 . 24 (sin  a)2[0.3365  + 0.3365 


cos  B cos  y + sin  B sin  y] 


(5.4) 
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The  Si-l^1  Si  DDF  is  related  to  the  distribution  V(a) 
of  the  Si-O-Si  bond  angle  a by  the  following  proportionality: 


ffr  I - XM 

G(rij3  T76'T~co~alT 


(5.5) 


n rl 

The  fourth  peak  (Si -2  0)  bond  length  involves  the  angle  3 

in  addition  to  the  angle  a.  Mozzi  and  Warren  found  that  a 
variation  of  3 such  that  3 takes  with  equal  probability  all 
values  from  0 to  360°  produces  a DDF  that,  yields  a PFD  con- 
sistent with  the  experimental  PFD  for  vitreous  silica. 

The  0-2nd  0,  and  Si-2nd  Si  distances  involve  the  angle  y 
in  addition  to  the  angles  a and  3.  Mozzi  and  Warren  found 
that  a variation  of  y similar  to  that  of  3 gave  the  best 
fit  to  the  experimental  PFD.  This  in  essence  means  that 
the  Si04  tetrahedra  in  vitreous  silica  can  assume  a random 
orientation  about  the  plane  containing  the  Si-0-Si  bond 
angle  a. 

The  x-ray  diffraction  patterns  of  the  heat-treated  sam- 
ples were  compared  with  patterns  of  mixtures  of  vitreous 
silica  powder  with  varying  percentages  of  powdered  cristo- 
balite.  The  areas  under  the  first  peaks  of  the  heat-treated 
sample  patterns  were  compared  with  the  corresponding  areas  of 
the  patterns  of  the  mixture  samples.  It  was  found  that  a 
close  match  in  the  areas  was  observable  between  the  19-hour 
sample  and  a mixture  containing  1-1/2  percent  cristobalite 
while  a match  was  obtained  between  the  30-1/2-hour  sample  and  a 


Figure  14.  Determination  of  the  percentage  crystallinity 
in  heat-treated  samples,  comparison  of  the 
first  diffraction  peak  areas. 

a)  Powder  heat-treated  19  hours  and  1-1/2 
percent  crystal  mixture. 


b)  Powder  sample  heat-treated  for  30-1/2 
hours  and  6 percent  crystal  mixture. 
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6 percent  mixture  sample  (Figs.  14  a and  b) . The  percentage 
of  crystallinity  in  the  heat-treated  samples  is  hence  small 
and  the  crystals  can  be  expected  to  be  well  distributed  in 
the  sample.  Hence  the  assumption  of  random  orientation 
could  be  made  for  obtaining  the  DDF  of  the  heat-treated 
samples  also. 

The  Si- 1st  Si  DDF  is  directly  connected  to  the  distri- 
bution V(cx)  of  the  Si-O-Si  bond  angle  a through  Eq.  (5.5). 
Hence  the  third  peak  DDF  could  be  calculated  directly  for 
various  assumed  values  of  the  V(ct)  distribution  or  vice 
versa.  The  Si-2n<^  0,  0-2nci  0 (to  be  called  the  5^  peak), 
and  Si-2nc^  Si  (henceforth  to  be  called  the  6t^1  peak)  bond 
length  calculations  involve  more  than  one  angular  distribu- 
tion. Hence  the  calculation  of  fourth,  fifth  and  sixth 
peak  DDF  is  more  involved.  The  DDF  of  peaks  4 to  6 were 
calculated  using  a computer  program.  The  procedure  for 
calculating  these  DDF  was  as  follows. 

Calculation  of  Fourth  to  Sixth  Peak  DDF 

Using  Eq.  (5.2),  the  maximum  and  minimum  values  of  the 
Si_2nd  q distance  were  determined  to  be  2.97  X and  4.27  A, 
respectively.  This  range  was  divided  into  intervals  of 
width  0.04  X.  A set  of  r; . values  was  calculated  by  vary- 
ing the  value  of  a between  120°  and  180°  in  5°  steps,  and 
for  each  value  of  a varying  8 from  0 to  180°  also  in  5° 
steps.  The  computer  program  then  classified  the  calculated 
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values  into  the  various  0.04  X intervals,  assigning  to 
each  r.  - value  the  weight  V(a)  of  the  corresponding  value 
of  a.  Considering  the  interval  of  r^  values  between,  say, 
2.97  X and  3.01  A,  there  may  be  , r-^  values  which  fall 
into  this  interval  wdien  the  value  of  a equals  r^ 

values  when  a equals  » and,  in  general,  Nn  r^j  values 
which  fall  into  this  interval  when  the  value  of  a is  equal 
to,  say,  ar . The  sum  of  the  products  NnV(an)  is  then  cal- 
culated for  each  of  these  intervals.  The  summation  Z N V (a  ) 

n 11  11 

increases  i\rith  r until  it  reaches  a maximum  value  in  a 
particular  interval.  During  the  course  of  the  present  in- 
vestigations, it  was  found  that  the  interval  in  which  the 
summation  attains  a maximum  value  shifted  from  4.17-4.21  X 
to  4.21-4.25  A as  a function  of  heat-treatment.  The  summa- 
tions Z N V(a  ) were  then  normalized  with  this  maximum  value 


n 


n v n 


as  one . 

The  fifth  and  sixth  peak  DDF  were  calculated  in  a 
similar  manner,  but  taking  into  account  the  variation  of 


the  third  angle  y for  each  of  the  set  of  values  of  a and  £. 
The  calculated  DDF  were  then  convoluted  with  the  correspond- 
ing pair  function,  Si-Si  in  the  case  of  the  third  peak,  Si-0 
for  the  fourth  peak,  0-0  in  the  case  of  the  fifth  peak  and 
so  on.  It  was  necessary  to  convolute  the  convoluted  pair 
functions  with  another  function  to  take  into  account  the 
contributions  of  thermal  vibrations  to  the  distribution  of 
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the  fourth,  fifth  and  sixth  peaks.  The  use  of  these  func- 
tions has  been  explained  by  Mozzi  and  Warren. 

Calculation  of  Theoretical  PFD 

The  pair  functions  which  have  been  convoluted  with  the 

proper  DDF  (and  thermal  vibration  functions  in  the  case  of 

the  fourth,  fifth  and  sixth  peaks)  were  then  substituted  in 

the  left-hand  side  of  Eq.  (3.2).  The  coordination  numbers 

N.  . were  the  same  as  determined  by  Mozzi  and  Warren.  Hence 
il 

the  value  of  the  theoretical  PFD  could  be  determined  as  a 
function  of  distance  from  a central  atom  taken  as  the  origin 
The  experimental  PFD  for  "T-D"  clear  fused  quartz  is  shown 
compared  with  the  PFD  calculated  by  assuming  a suitable  set 
of  DDF  (Figs.  15a  and  15b).  The  individual  calculated  peaks 
for  "T-D"  clear  fused  quartz  are  shown  in  Fig.  15a,  while 
Fig,  15b  shows  the  composite  curve  obtained  by  adding  to- 
gether the  individual  contributions  through  the  range  0 to 

O 

7.1  A.  It  was  found  that  the  third  peak  DDF  were  similar 
to  that  obtained  by  Mozzi  and  Warren.  The  comparison  of 
the  experimental  and  the  calculated  (total)  PFD  for  the 
heat-treated  samples  is  shown  in  Figs.  16  and  17.  It  is 
clear  that  a fairly  good  fit  is  obtained  for  the  first 
three  peaks  of  the  PFD,  while  the  fit  is  not  so  good  from 
the  fourth  peak  on.  The  goodness  of  fit  was  judged  by  the 
following  criteria: 


Figure  15a.  Comparison  of  the  experimental  PFD  of 
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Figure  15b.  Comparison  of  the  experimental  PFD  of  "T-D"  clear 
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1.  the  peak  positions  of  the  experimental  and 
calculated  theoretical  PFD  should  coincide, 
and 

2.  the  DDF  were  altered  and  various  PFD  calcu- 
lated for  any  given  set  of  Sfflax  and  a values 
(corresponding  to  the  experimental  values) . 

The  set  of  DDF  which  gave  the  closest  match 
out  to  the  third  peak  of  the  experimental  PFD 
was  taken  to  represent  the  structure  under 
consideration . 

The  lack  of  fit  from  the  fourth  peak  on  is  due  to  the  fact 
that  contributions  from  peaks  beyond  the  sixth  peak  have 
not  been  taken  into  account.  Mozzi  and  Warren  have  shown 

„ J O 

that  a Si -3  0 peak  should  exist  at  about  6.4  A.  The  con 

tribution  of  this  peak  as  well  as  subsequent  peaks  is  neg- 
lected in  the  PFD  calculations.  However,  it  can  be  seen 
that  the  positions  and  the  shapes  of  the  fourth  and  the 
combination  of  the  fifth  and  sixth  peaks  in  the  calculated 
PFD  are  about  right. 

Heat -Treatment.  Induced  Changes  in  the 
Vitreous  Silica  DDF 

The  corresponding  V(a)  and  DDF  curves  are  shown  in 
Figs.  18  to  21.  A theoretical  DDF  for  a sample  heat- 
treated  for  a time  longer  than  30-1/2  hours  is  also  in- 
cluded in  Figs.  18  to  21  for  comparison.  An  examination 
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Figure  19.  Variation  of  Si-2  0 DDF  with  heat -treatment 
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Figure  20.  Variation  of  0-2  0 DDF  with  he at -treatment . 
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Figure  21.  Variation  of  Si-2  1 “ Si  DDF  with  heat -treatment . 
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of  Figs.  18  through  21  reveals  pronounced  changes  occurring 
in  the  V(a)  and  Si-lst  Si  DDF  as  a function  of  heat-treat- 
ment and  less  dramatic  changes  occurring  in  the  DDF  of  the 
fourth,  fifth  and  sixth  peaks.  The  changes  in  the  DDF  half- 
breadths  are  summarized  in  Table  6 (which  also  gives  the 
total  breadths  of  the  DDF)  and  Fig.  22.  AX  in  Fig.  22  repre- 
sents the  breadth  at  half -maximum  of  the  V(a)  as  well  as  the 
corresponding  DDF  curves  for  the  heat-treated  samples,  while 
AX (glass)  gives  the  corresponding  parameters  for  "T-D"  clear 
fused  quartz.  The  ratio  AX/AX (glass)  for  the  crystalline 
sample  is  represented  as  a single  dashed  line  near  zero. 

This  ratio  is  represented  to  be  slightly  greater  than  zero 
as  the  DDF  would  still  have  a slight  breadth  due  to  temper- 
ature vibrations  and/or  imperfections  in  the  crystals  and 
so  on.  It  is  clear  that  the  Si-lst  Si  DDF  is  approaching 
the  crystalline  value  faster  than  the  other  DDF. 

Interpr etation  of  Heat -Treatment  Results 

The  bulk  heat- treatment  results  show  that  there  is  no 
noticeable  relaxation  in  the  structure  of  vitreous  silica 
as  a function  of  heat- treatments . The  DDF  of  heat-treated 
powder  samples,  on  the  other  hand,  show  changes  as  a result 
of  heat-treatment.  The  changes  could  be  analyzed  in  the 
foil o w i n g manner. 

Interpretations  of  the  powder  heat -treatment  results 
as  bulk  rather  than  localized  structural  alterations  require 
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Table  6 

Breadth  of  the  DDF  for  Various  Ion  Distances 
in  Unheat-Treated  Vitreous  Silica 


Si- 1st  Si  Si - 2nd  0 0- 2nd  0 Si-2nd  Si 


Half- 

breadth  0.32  A 0.4  A 0.7  1 1.2  A 


Total  0 

breadth  0.44  A 1.4  A 3.1  A 3.5  1 
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Figure  22.  Changes  in  the  half-maximum  breadths  of  SiO?  DDF  with 
heat  - treatment . 
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that  surface  crystallization  of  the  powder  does  not.  pro- 
duce comparable  changes  in  the  DDF.  Transmission  electron 
microscopy  of  the  heat-treated  powders  reveal  the  presence 
of  crystallites  in  the  surface  of  the  particles  (Figs.  23 
and  24)  . The  presence  of  a crystallite  on  the  surface  of 
an  amorphous  particle  is  shown  in  the  electron  diffraction 
pattern  o£  a particle  from  the  19-hour  sample  (Fig.  23). 

A particle  from  the  30 -1/2 -hour  sample  appears  to  be  fully 
crystallized.  We  have  seen  that  the  analysis  of  the  x-ray 
Bragg  peaks  of  the  heat-treated  powder  samples  reveals  the 
presence  of  1-1/2  percent  and  6 percent  crystals  in  the  19- 
and  30-1/2-hour  heat-treated  samples,  respectively  (Figs. 

14  a and  b)  . Synthetic  mixtures  containing  1-1/2  percent 
and  6 percent  crystals  were  analyzed  by  the  PFD-DDF  method. 
The  PFD  of  the  results  obtained  are  shown  in  Figs.  25-27. 

The  PFD  of  the  powder  heat-treated  for  19  hours  at  1200°C  is 
shown  compared  with  the  PFD  of  the  corresponding  synthetic 

sample  in  Fig.  25.  The  S values  for  the  19 -hour  and  the 
r max 

synthetic  samples  are  1.85  A \ and  1.80  A \ respectively. 
This  factor  should  cause  the  19 -hour  PFD  peaks  to  be  sharper 
than  the  synthetic  PFD  peaks,  even  if  the  structure  of  the 
heat-treated  and  the  synthetic  samples  were  the  same.  For 
example,  the  DDF  obtained  for  the  powder  sample  heat-treated 
for  19  hours  were  then  used  to  calculate  the  PFD  using  the 
S value  appropriate  to  the  1-1/2  percent  crystal  samples. 
The  calculated  PFD  is  compared  to  the  experimental  PFD  in 


83 


A 

+j 

s 

o 


rt 

K) 

CVI 

<L> 

3 

W> 

•H 


sample  heat-treated  19  hours  at  1200°C  (magnification  50,000X). 
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Figure  23b 


Electron  diffraction  pattern  of  particle 
in  Fig.  23a. 
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Figure  24b.  Electron  diffraction  pattern  of  particle 
in  Fig.  24a. 


£' 


ID 

P 

P 

O 

rP 

C7i 

rH 

<D  ' 
4.1  • 

CtS  (D 
0)  P 
P 3 
■p  p 
• X 
P -H 
nS  £3 
o 

,P  .H 
cC 

<u  P 
rp  1/1 

p.  >-> 

e p 

as  O 


P P 

O <D 
r<3  O 

> i_i 


o c 

p,  p, 


<4-1  (SI 


t-J.  1 
Cl—  rH 

Pi 

•P 
tH  P 
rt  d 
P 

P u 
0)  o 
£ O 
•H  O 
P (s] 
4)  t — 1 

Ph 
X P 
yo  rt 


LO 

(SI 

P 

P 

P 

60 

•H 

IP 


19  HOURS  { S„AX~  !-85 


88 


09  ZL  V9  09  8’  > C>  Z£  i-2  9!  90  00 

1 * i 1 1 < J i_ f-00 


igure  26.  Comparison  of  the  experimental  PFD  for  the  1-1/2  percent 
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Fig.  26.  The  fit  is  fairly  good,  except  for  a shift  in  the 
position  of  the  calculated  third  peak.  Differences,  if  any, 
between  the  19-hour  heat-treated  and  mixture  samples  are  too 
small  to  be  resolved  by  the  PFD-DDF  approach. 

The  PFD  of  the  30-1/2-hour  heat-treated  sample  fits 
rather  closely  with  the  corresponding  synthetic  PFD  (Fig. 
27).  The  DBF  analysis  confirms  this  fit.  The  observations 
indicate  that  there  may  be  a small  volume  relaxation  effect 
prior  to  the  onset  of  crystallization  in  vitreous  silica, 
but,  if  so,  it  is  barely  within  the  detection  limits  of  this 
method.  The  observed  slight  changes  may  also  be  explainable 
as  due  to  the  local  changes  in  the  structure  of  vitreous 
silica  at  the  surfaces  in  contact  with  the  crystallites. 
Because  of  changes  in  the  density, as  the  crystallites  grow 
they  fracture  from  the  surface  and  hence  cannot  influence 
the  surface  structure  of  the  remaining  particle  in  a manner 
similar  to  the  initial  stages  of  crystallization. 


Effect  of  Water  Content,  on  the 
~ Structure  of  Vitreous  Silica 


The  effect  of  water  content  on  the  prop 
silica  has  been  briefly  mentioned  in  Chapter 
ing  the  marked  variations  in  properties  such 
[Hetherington  et  al . (1964 ) ] , crystallization 
[Wagstaff  et  al . (1964)],  optical  absorption 


erties  of  fused 
II.  Consider- 
as  viscosity 
kinetics 
[Stolen  et  al. 
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(1971)],  and  dielectric  relaxation  [Owen  and  Douglas  (1959)] 
with  variation  in  water  content,  it  may  be  expected  that 
water  should  have  a noticeable  effect  on  the  structure  of 
vitreous  silica.  It  is  generally  agreed  that  water  is 
present  in  combined  form  in  fused  silica.  Each  w ater  mole- 
cule is  incorporated  in  the  silica  structure  by  the  fracture 
of  a Si-Q-Si  bridge  and  the  formation  of  two  silanol  groups, 
a reaction  which  can  be  represented  as  follows  [Adams  and 
Douglas  (1959),  Rawson  (1967)]. 


\ 

— — Si 0 — 


H.OH 


If  this  mechanism  for  the  incorporation  of  water  in 
vitreous  silica  is  true,  then  this  should  cause  rupture  of 
Si-O-Si  bonds  and  thus  should  increase  the  Si-lst  Si  dis- 
tances involved.  Hence  the  effect  of  water  content  should 
be  most  pronounced  on  the  third  peak  DDF  and  should  be 
noticeable  with  increasing  water  content  (see  also  Appendix 
4)  . 

PEP  Analysis  of  Samples  with  Varying 
Water  C o n tent s 

La  the  present  investigation,  fused  silica  samples 
covering  a range  of  water  contents  over  three  orders  of 
magnitude  were  studied.  IR  Vitreosil  contains  0.0003  percent, 
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V.itreosil  0()'5,  0.04  percent,  while  "T-D"  clear  fused  quartz 
and  Spectrofnil-B  both  contain  0.12  percent  (all  weight  per- 
cent as  OH) . It  should  be  noted  that  all  the  samples  except 
Spectrosil-B  are  produced  by  fusion  of  SiO?  rock  crystal. 

The  scattered  intensity  in  absolute  units  per  unit  of  compo- 
sition for. .-.these  four  samples  is  shown  in  Fig.  28.  The  in- 
* 

tensity  curves  of  all  these  samples  except  for  Spectrosil-B 

are  nearlyiidentical . The  Spectrosil-B  curve  shows  a 

sharper  first  maximum  and  more  oscillations  at  high  S values 

than  the  others.  The  corresponding  PFD  are  compared  in  Fig. 

29.  The  S for  the  "T-D"  clear  fused  quartz  was  forced  to 

coincide  with  the  S of  the  other  samples  to  make  direct 

max  1 

comparison  possible.  The  PFD  of  all  the  samples  are  nearly 
identical  except  that  the  Spectrosil-B  PFD  has  a.  narrower 
third  peak  and  the  fourth  and  fifth  peaks  (containing  both 
the  0-2nc^  0 and  Si-2nc^  Si  interactions)  are  sharper  and 
narrower  for  the  Spectrosil-B  PVD  than  the  others.  A DDF 
analysis  of  the  Spectrosil-B  PFD  was  performed  using  the  DDF 
for  the  powder  heat-treated  for  19  hours.  A fairly  good  fit 
was  observed  between  the  experimental  and  calculated  PFD, 


thus  indicating  the  possible  presence  of  a crystalline 
structure  in  the  Spectrosil-B  sample.  The  following  con- 
clusions can  be  drawn  from  these  results. 


T D FUSED  QUARTZ 
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Figure  28.  Normalized  intensity  curves  for  samples  containing  various  i 


4-> 

£ 

a; 

o 

E 

0 

P< 

W 

3 

O 

•H 

E 

cti 

> 

be 
£ 
•H 
— ( 
*-h 

•H 

rt 

4-J 

£ 

o 

O 


to 

0 

rH 

£•< 

£ 

rt 

10 

<p 

o 

n 

PH 

pH  * 

E 

rH  O 

Ct!  P 
P E 
E 
0 

£ m 

•H  O 
E 

0 V) 
pH  0 

rt  bC 
W 3 


CTl 

CM 

0 

E 

3 

bo 

•rH 

PH 


T 0 FUSED  QUARTZ 


98 


q 

co 


OJ 


u> 


. q 

SO 


CO 

< 


_ « 

80 


Ol 


<0 


! ° 


o 

d 


o 

0 

•o 

1 


ccj[ 


99 


Interpretation  of  the  PFD  Results 

Vitreous  silica  samples  produced  by  fusion  of  rock  crys- 
tal have  a similar  structure  although  their  water  contents 
vary  from  0.0003  percent  to  0.12  percent.  Spectrosil-B,  on 
the  other  hand,  behaves  as  though  it  were  more  crystalline 
than  the  other  samples.  Spectrosil-B  is  the  product  of 
vapor-phase  hydrolysis  of  pure  SiCl^  in  a flame.  Hence  the 
existence  of  "residual  crystallinity"  should  be  less  probable 
in  the  case  of  Spectrosil-B  than  for  the  samples  produced  by 
fusion  of  quartz  crystal.  However,  the  study  of  Owen  and 
Douglas  (1959)  on  the  electrical  properties  of  the  various 
Thermal  Syndicate  fused  silica  samples  should  be  considered 
here.  Owen  and  Douglas  were  studying  the  dielectric  losses 
in  the  various  grades  of  silica  as  a function  of  temperature 
near  the  quartz  inversion  temperature  (573°C).  They  found 
that  the  dielectric  loss  curves  of  all  the  samples  tend  to 
have  a maximum  at  about  575°C.  This  maximum  was  most  pro- 
nounced in  the  case  of  Spectrosil-B.  Owen  and  Douglas  con- 
cluded that  this  was  evidence  for  a quartz-like  transition 
in  fused  silica.  It  is  surprising  that  this  should  be  most 
pronounced  in  the  case  of  Spectrosil-B,  which,  so  far  as  is 
known,  has  never  passed  tnrough  a quartz  stage.  It  is  known 
[Kracek  (1951)]  that  water  can  act  as  a mineralizer  in  fused 
silica.  It  is  possible  that  the  high  (0.12  weight  percent) 
water  content  in  Spectrosil-B  could  have  promoted  formation 
of  crystalline  regions  (presumably  quartz).  Hetherington 


100 


et  al.  (1964)  have  noted  that  the  water  is  more  uniformly 
distributed  in  a synthetic  vitreous  silica  (Spectrosil -8) 
compared  to  the  flame -fused  grades  where  a gradient  was 
observed  in  each  grain.  This  may  explain  why  the  same  water 
content  in  the  flame-fused  "T-D"  clear  fused  quartz  does  not 
have  such  a widespread  and  pronounced  "mineralizing"  effect 
as  in  Spect rosil-B . 


Effect  of  Neutron  Irradiation  on 
the  Structure  of  Vitreous  Silica 

The  normalized  unmodified  intensity  data  for  the  Opto- 

sil  sample  before  and  after  a fast  neutron  irradiation  to  a 

7 0 2 

total  flux  of  about  10'  /cm  are  shown  in  Fig.  30.  The 
first  peak  of  the  diffraction  pattern  is  broader  and  con- 
siderably reduced  in  intensity  as  a result  of  the  irradia- 
tion. The  peaks  at  S values  of  0.7  X ^ and  0.79  X ^ merge 
into  a single  peak  at  about  0.77  X ^ . The  other  features 
of  the  Intensity  curves  are  essentially  unaltered  as  a re- 
sult of  Irradiation. 

Chan ges  in  the  PFD  Caused  by 
Neutron  Irradiation 

The  corresponding  experimental  PFD  are  compared 
in  Fig.  31.  A study  of  Fig.  31  shows  that  heavy  irradiation 
has  resulted  in  the  following  changes: 

1)  a shift  in  the  position  of  the  first  peak  from 
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Figure  30.  Normalized  intensity  curves  for  Optosil  samples  before  and  after 
heavy  neutron  irradiation. 
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1.62  A to  about  1.65  A. 

2)  a shift  in  the  position  of  the  second  peak  with 
a merging  of  the  second  and  third  peaks.  The 
shoulder  between  the  second  and  third  peaks  is 
no  longer  observable. 

3)  the  fourth  peak  seems  to  be  affected  little  if 
at  all,  while 

4)  the  fifth  peak  has  gone  down  in  value.  This  is 
in  spite  of  the  fact  that  the  density  and  there- 
fore the  average  electron  density  used  to 
evaluate  the  experimental  PFD  is  about  3 percent 
higher  for  the  irradiated  sample. 

DDF  Analysis  of  the  Irradiated  Sample 

A detailed  DDF  analysis  was  performed  on  the  experi- 
mental PFD  of  the  two  samples.  The  DDF  of  the  Optosil  sam- 
ple before  irradiation  is  essentially  the  same  as  the  DDF 
of  the  !fT-D"  clear  fused  quartz  sample.  The  DDF  of  the 
irradiated  sample  has  the  following  characteristics: 

1)  the  first  peak  DDF  is  shifted  to  slightly  higher 
values  of  r with  the  maximum  at  1.65  A instead 
of  1.62  A. 

2)  the  second  peak  DDF  had  also  to  be  shifted  to 
higher  values  with  the  distribution  of  atoms 
increasing  in  the  region  from  2.65  to  2.9  A 
instead  of  being  symmetrical  about  2.65  A and 
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falling  off  to  zero  at  2.85  X,  as  in  the  case 
of  the  corresponding  DDF  before  irradiation. 

3)  pronounced  changes  are  observed  in  the  third 
peak  DDF.  It  is  broader  than  the  corresponding 
DDF  before  irradiation,  and  extends  over  a wider 
range  of  r values  (Fig.  32  a). 

The  observed  changes  are  explainable  if  the  assumption 
is  made  that  the  first  peak  is  shifted  from  1.62  X to 

O 

1.65  A.  It  would  be  expected  that  if  the  tetrahedral  coordi- 
nation is  unaffected  the  second  peak  should  now  occur  at 
about  2.70  X instead  of  at  2.65  X.  This  seems  to  be  con- 
firmed by  observation  2 above.  Knowing  the  third  peak 
(Si-lst  Si)  DDF,  G(r^),  a and  the  Si-lst  0 distance,  the 
distribution  V(a)  of  the  5i-0-Si  bond  angle  a in  the  irradi- 
ated sample  can  be  calculated  from  the  proportionality 

V (a)  or  G (ri j ) 1.65  cos  a/2  (5.6) 

where  a is  given  by  the  relation 

sin  a/2  = r^/3.30  (5.7) 

The  V(a)  calculated  from  the  experimental  DDF  using  Eqs. 

(5.6)  and  (5.7)  is  shown  compared  to  the  V(a)  curve  before 
irradiation  (Fig.  32  b) . The  V(a)  curve  has  shifted  to 
lower  values  of  a and  the  maximum  in  the  V(a)  curve  occurs 
at  about  138°  (instead  of  the  144°  observed  for  the  un- 
irradiated sample) . 
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Using  the  increased  value  of  the  Si-0  distance,  a new 
set  of  DDF  for  the  peaks  4 to  6 was  calculated.  These  DDF 
were  a set  used  to  calculate  the  theoretical  PFD  for  the 
irradiated  sample.  The  experimental  and  calculated  .PFD  are 
shown  compared  in  Fig.  33.  It  is  clear  that  the  assumed  DDF 
values  give  rise  to  a PFD  which  compares  well  with  the  ex- 
perimental PFD.  All  the  other  assumptions  which  were  tried 
failed  to  give  such  a close  fit.  The  DDF  of  peaks  4 to  6 
are  shown  in  Figs.  34  to  36.  A general  widening  of  the 
•atomic  distribution  as  a result  of  the  irradiation  treat- 
ment is  noticeable. 


The  Structure  of  Neutron- Irradiated 
Vitreous  Silica 

The  PFD-DDF  results  indicate  that: 

s t 

1)  there  is  an  increase  in  the  mean  Si-1'  0 

distance  from  1.62  A to  1.65  A. 

2)  the  tetrahedral  units  of  vitreous  silica  are 
essentially  preserved  even  after  heavy  neutron 
irradiation . 

s t 

3)  the  Si-1  Si  distances  tend  to  be  distributed 
over  a wider  range  of  r values  with  a concomi- 
tant shift  in  the  distribution  of  the  Si-O-Si 
bond  angles  towards  lower  a values. 


5 000 


108 


Figure  33.  Comparison  of  the  experimental  PFD  for  the  irradiated  sample  with 
the  calculated  PFD  (composite) . 
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Figure  34.  Variation  of  the  Si-2  0 DDF  as : a result  of  neutron  irradiation. 
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Figure  35.  Variation  of  the  0-2  0 DDF  as  a result  of  neutron  irradiation. 
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4)  the  other  DDF  also  tend  to  be  distributed  over 
wider  ranges  of  r values  as  a result  of  irradi- 
ation. 

Primak  (1958)  has  discussed  the  effect  of  heavy  neutron 
irradiation  and  subsequent  annealing  of  vitreous  silica  on 
properties  such  as  density  and  refractive  index.  His  con- 
clusions are  instructive.  He  proposes  that  a thermal  spike 
mechanism  is  responsible  for  the  observed  changes  on  irradi- 
ation of  vitreous  silica.  Simon  (1960)  and  Bale,  Shepler 
and  Gibbs  (1970)  also  propose  that  a thermal  spike  mechanism 
is  responsible  for  the  compaction  of  and  change  in  the  struc- 
ture of  vitreous  silica.  The  temperature  in  the  interior  of 
the  thermal  spikes  is  expected  to  rise  as  high  as  2500°C. 
However,  this  region  will  be  essentially  quenched,  freezing 
the  high  temperature  structure  within  times  of  the  order  of 
10  ^ sec  [Primalc  (1958)].  The  number  and  distribution  of 
the  thermal  spikes  may  be  expected  to  be  a maximum  at  satu- 
ration irradiation.  It  would  also  be  expected  that  the 
deformations  of  the  tetrahedra  themselves  will  be  slight 
[Primak  (1958)],  but  greater  deformations  of  these  configu- 
rations writh  respect  to  one  another  can  occur  due  to  the 
collision  processes  occurring  during  the  neutron  bombard- 
ment. This  seems  to  occur  in  such  a way  as  to  make  the 
structure  more  compact  with  respect  to  the  SiO^  tetrahedra 
to  compensate  for  the  increased  Si-1'  0 distances.  This 
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compacted  structure  must  then  be  frozen  in  as  the  thermal 
spike  chills. 

Lorch  (1969)  has  also  observed  an  increase  in  the  mean 
interatomic  distances  in  his  sample  of  Spectrosil  as  a re- 
sult of  neutron  irradiation.  The  PFD-DDF  analysis  is  much 
more  powerful  and  the  present  results  show  that  a wide  dis- 
tribution of  interatomic  distances  can  be  expected  as  a 
result  of  neutron  irradiation.  Lorch’ s results  seem  to  be 
in  qualitative  agreement  with  the  present  PFD-DDF  results 
of  the  structure  of  irradiated  vitreous  silica. 


CHAPTER  VI 


SUMMARY  AND  CONCLUSIONS 

Mozzi  and  Warren’s  new  PFD-DDF  analysis  provides  a 
basis  for  a quantitative  description  of  the  structure  of 
vitreous  silica.  The  DDF  of  the  various  interatomic  dis- 
tances in  vitreous  silica  provide  such  a quantitative  de- 
scription. This  approach  is  more  precise  than  the  early  RDF 
analysis  methods.  The  DDF  by  themselves,  or  a characteris- 
tic feature  of  the  DDF  such  as  the  breadth  at  half-maximum 
of  the  DDF  curves,  can  be  used  to  describe  a given  struc- 
tural state  of  silica.  The  present  investigation  has  used 
the  DDF  and  their  breadths  (at  half-maxima)  to  give  a pre- 
cise quantitative  description  of  the  structural  changes  in 
vitreous  silica  as  a function  of  heat -treatment , water  con- 
tent and  heavy  neutron  irradiation. 

The  following  conclusions  can  be  drawn  from  the  results 
of  the  present  investigation. 

1.  Heat-treatment  of  bulk  fused  silica  at  temperatures 
between  1400°C  and  1661°C  does  not  cause  a noticeable  relax- 
ation in  the  structure. 

2.  Heat- treatment  of  powdered  silica: 

a)  a definite  time  period  is  necessary  before 
crystallization  in  powdered  silica  can  be 
observed . 
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b)  the  PFD  of  the  powder  sample  heat-treated 
for  19  hours  at  1200°C  showed  a slight  dif- 
ference from  the  PFD  of  the  corresponding 
crystalline  mixture.  The  DDF  analysis  re- 
veals that  differences,  if  any,  in  the  struc- 
ture of  the  crystal  line  mixture  and  the 
heat-treated  powder  sample  are  just  at  the 
limit  of  resolution  of  the  present  tech- 
niques . 

c)  the  PFD-DDF  of  the  powder  heat-treated  for 
30-1/2  hours  at  1200°C  is  essentially  similar 
to  the  PFD-DDF  of  the  corresponding  crystal- 
line mixture. 

d)  from  (b)  and  (c)  above  it  would  seem  that 
partially  crystallized  silica  powder  is  simi- 
lar to  a mixture  of  cristobalite  and  vitreous 
silica. 

e)  a detailed  DDF  analysis  has  been  performed 

on  the  heat-treated  powder  samples,  and  shows 

s t 

considerable  differences  in  the  Si-1  Si  and 
Si-Z1105  Si  DDF  as  a result  of  heat -treatment , 
while  less  pronounced  changes  occur  in  the 
other  DDF. 

3.  The  PFD  of  "T-D"  clear  fused  quartz,  Vitreosil  066 
and  IR  Vitreosil  are  closely  similar  to  one  another,  while 
the  PFD  of  Spectrosil-B  indicates  a more  ordered  structure. 
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The  first  three  specimens  span  a range  of  water  contents 
from  0.0003  to  0.12  weight  percent,  and  are  all  produced  by- 
fusion  of  rock  crystal.  Spectrosil-B,  on  the  other  hand, 
is  a synthetic  fused  silica  with  the  same  water  content  as 
"T-D"  clear  fused  quartz.  These  results  are  in  qualitative 
agreement  with  the  dielectric  relaxation  results  of  Owen  and 
Douglas  (1959) . 

4.  The  PFD-DDF  analysis  of  an  Optosil  sample  after 
heavy  neutron  irradiation  reveals  an  increase  in  the  dis- 
tribution of  interatomic  distances  in  fused  silica  as  a 
result  of  irradiation,  with  a shift  in  the  Si-O-Si  bond 
angles  towards  lower  values  of  a.  The  results  of  the  pres- 
ent investigation  confirm  the  main  details  of  the  mechanism 
for  radiation  compaction  of  vitreous  silica  proposed  by 
Primak  (1958)  . 

5.  Silicic  acid  containing  88  percent  SiC^  and  the 
remainder  essentially  water  was  also  investigated.  The  PFD 
of  this  silicic  acid  shows  a considerable  change  compared 
with  the  PFD  of  "T-D”  clear  fused  quartz.  The  differences 
are  due  to  the  presence  of  water  in  the  structure  as  com- 
bined water  which  disrupts  the  Si04  network,  changing  the 
nearest  neighbor  distance  distributions  and  coordinations. 


CHAPTER  VII 


SUGGESTIONS  FOR  FUTURE  INVESTIGATORS 

1.  During  the  course  of  the  present  investigation,  the 
bulk  heat -treatments  were  carried  out  at  1400°C,  1450°C, 

1 595 °C  and  I661°C  only.  The  heat -treatments  should  be  ex- 
tended to  other  temperatures  starting  with  1200cC  to  deter- 
mine if  any  structural  relaxation  can  occur  in  vitreous 
silica  at  temperatures  closer  to  the  silica  transformation 
range . 

2.  Seeding  with  colloidal  metal  particles  such  as  Pt 
dissolved  in  Si02  melt  should  be  tried.  The  precipitation 
of  the  colloidal  particles  should  help  in  a more  uniform 
nucleation  in  the  bulk  of  silica.  Structural  changes,  if 
any  should  occur  on  nucleation  and  crystallization  in  the 
bulk  of  silica,  will  be  more  widespread  and  capable  of  be- 
ing followed  successfully  by  the  PFD-DDF  approach. 

3.  The  DDF  approach  should  be  extended  to  binary  sili- 
cate glasses  containing  network  modifying  cations.  It  might 
be  expected  that  the  x-ray  diffraction  pattern,  and  hence 
the  structure,  is  a function  of  heat -treatment  in  such  com- 
positions and  should  be  amenable  to  investigations  using 
the  present  approach. 
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4.  X-ray  small  angle  scattering  studies  should  be  made 
concurrently  with  the  large  angle  studies  to  detect  the 
presence  of  nuclei  in  the  early  stages  of  crystallization 
(especially  for  2 and  3 above) . 

5.  Transmission  electron  microscopy  should  also  be 
used  to  follow  the  initial  stages  of  the  crystallization  of 
vitreous  silica.  Zarzycki  (1970)  has  used  a hot  stage  elec- 
tron microscope  to  study  changes  in  filaments  of  vitreous 
silica  due  to  heating  in  the  electron  beam  and  has  observed 
formation  of  domains  in  the  filaments.  It  would  be  of  inter- 
est to  see  if  such  domains  are  formed  in  the  bulk  samples, 
and,  if  so,  at  what  level  would  these  affect  the  x-ray  dif- 
fraction results.  This  would  especially  be  a fruitful 
analysis  in  the  case  of  vitreous  silica  samples  seeded  with 
colloidal  particles  to  ensure  widespread  nucleation  in  the 
bulk  of  the  sample. 

6.  The  use  of  powder  should  provide  a viable  method 
for  following  surface  structural  changes  as  a result  of  cor- 
rosion, ion-exchange,  and  other  surface  treatments  which 
have  a definite  effect  on  the  properties  of  glasses.  It 
might  be  expected  that  the  use  of  fine  powder  would  provide 
a large  enough  surface  area  so  that  structural  changes,  if 
any,  can  be  easily  determined  by  the  PFD-DDF  approach. 

7.  The  structural  changes  occurring  in  vitreous  silica 
as  a result  of  heavy  neutron  irradiation  have  been  analyzed 
in  the  present  investigation.  It  would  also  be  instructive 
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to  follow  the  structural  changes  using  the  PFD-DDF  approach 
as  a function  of  varying  levels  of  neutron  irradiation. 

8.  The  structural  changes  occurring  as  a result  of 
annealing  the  neutron  irradiated  samples  should  also  be  fol- 
lowed using  the  PFD-DDF  techniques.  This  should  clarify  the 
structural  aspects  of  neutron  irradiation  damage  in  vitreous 
silica.  This  analysis  could  be  extended  to  glasses  other 
than  vitreous  silica  as  well. 

9.  Improvements  should  be  made  to  the  x-ray  diffrac- 
tion system  used  in  the  present  investigation.  The  use  of 
Mo  radiation  in  the  present  investigation  limited  the  range 
of  S values  in  which  accurate  intensity  data  were  available . 

An  experimental  set-up  using  rhodium  or  silver  radiation  with 
primary  beam  monochromator  and  the  Warren-Mavel  fluorescence 
technique  to  experimentally  eliminate  the  Compton  scatter- 
ing should,  be  tried  to  give  accurate  intensity  data  out  to 
higher  values  of  S.  This  would  enhance  the  resolution  pos 
sible  in  the  PFD  analysis  method  and  make  the  DDF  analysis 
more  precise. 

10.  In  cases  where  it  is  not  possible  to  experimentally 
eliminate  Compton,  the  Compton  scattering  should  be  experi- 
mentally determined  as  a function  of  S.  This  is  necessary 
as  the  tabulated  values  of  Compton  are  often  suspect,  and 
the  Compton  scattering  is  slightly  sensitive  to  the  type  of 


bonding  in  the  glass,  etc. 
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11.  Infra-red  studies  should  be  made  of  the  structure 
of  glass  as  a function  of  the  various  variables  mentioned 
above  and  the  DDF  and  infra-red  results  combined  to  give  a 
better  picture  of  the  structure. 

12.  The  structure  of  boric  oxide  and  boro-silicate 
glasses  should  also  be  investigated  using  the  approach  out- 
lined in  the  present  work.  The  question  of  degree  of  ion- 
icity  in  the  boric  oxide  case  is  interesting  and  should  be 
attempted . 


APPENDICES 


APPENDIX  1 


NORMALIZATION  OF  COPPER  INTENSITY  DATA 

Let  Ij,  (in  counts  per  second)  denote  the  measured  in- 
tensity, A the  air-scattering  contribution,  B the  background 
noise,  and.  p.f.  the  polarization  factor  for  primary  beam 
monochromated  Cu.  We  want  to  find  a constant  K such  that 


K[ 


p.f. 


I 

e .u. 


N 


+ Z 

u . c . 


i . 
J 


(A-l.l) 


where  I /N  represents  the  intensity  in  electron  units 
scattered  by  one  unit  of  composition,  in  this  case  SiC^, 
omitting  the  small  angle  scattering  contribution,  and  i ^ is 
the  Compton  scattering  from  atom  j . 


Now , 


i(10 


I 

e . u. 


N 


Z f 4 2/g2  (k) 

u.  c . J 


[ from  Eq . (3.5)], 


IM-(A+B)  . 2 

hence,  i ($0  = K[ — — — ? ] - Z i-  - I f-  (A- 1.2) 

P* r ■ u.c.  J u.c.  J 


We  know  from  Warren  (1969b)  that  for  vitreous  silica 


k 

/ 


i(k) 


2.2 

.2  -a  k 
k e 

g2(k) 


dk 


- 2 TT  p E 
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Z.  + 29 
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(A-1.3) 
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Substituting  for  i(k)  from  Eq.  (A- 1.2)  we  get 


k ,2  -a2k2  I..-CA+B)  k ,2  -a2k2 

irrITlkG  f M l r PI  k 6 Tr  C -P  ^ j-  A ^ 1 

k f — 7 L — zr-p Jak  - / — J [ E (£•  -1-JJ 

g (k)  p-x‘  ° g2  fk)  u.c.  J 3 


- 2tt  p E Z.  + 29  (A- 1.4) 

u.c . 3 

IM" CA+B) 

Experimental  values  of  — — j were  available  only  down  to 

a 20  value  of  10°  (S  ~ 0.11  A"').  A smooth  extrapolation 
of  the  experimental  values  could  be  performed  to  an  S value 
of  zero  [Warren  (1969b)]. 

A computer  program  was  written  to  calculate  the  inte- 

0 -1 

grals  involved  in  Eq.  (A-1.4).  "k  " was  taken  as  7.5  A 
(Sm  = 1.2  X~^) . Knowing  p£  and  Z.,  the  normalization  con- 
stant K was  computed  and  the  normalized  unmodified  intensity 
in  electron  units  per  unit  of  composition  determined  using 
Eq.  (A- 1.1) . The  program  gives  the  following  output: 

1.  printout  of  the  values  of  the  unmodified  in- 
tensity in  absolute  units  (electron  units  per 
unit  of  composition)  as  a function  of  S in 
steps  of  0.012  X \ 

2.  punched  card  output  in  a format  suitable  for 
use  in  the  subsequent  program  for  calculating 
the  PFD,  and 

a plot  of  intensity  as  a function  of  S from 

0 to  1.2  X"1. 


3. 


APPENDIX  2 


KAPLOW  CORRECTION  PROCEDURE 


Let  I (k)  be  the  true  scattered  intensity  in  electron 
units  per  unit  of  composition  given  as  a function  of  k 
(2irS)  . A function  F(k)  can  be  defined  such  that 


F(k)  = k[- 


I(K) 


I 

u . c , 


- 1]  = /°°  G(r)  sin  kr  dr 
l o 


(A-2.1) 


and 


G(r)  = 2/'tt  F(k)  sin  rk  dk  (A-2.2) 

where  G(r)  = 4ivr [ p (r) - p ] (A-2.3) 


pQ  is  the  average  atomic  density  over  the  sample  volume, 
and  p(r)  is  the  average  atomic  density  distribution  about 

O 

an  arbitrarily  chosen  origin  to  r values  of  about  10  A.  In 
ether  words,  G(r)  is  one  way  of  representing  an  RDF. 

The  experimental  intensity  values  and  hence  the  result- 
ing PFD(RDF)  are  subject  to  the  following  errors: 

1.  experimental  errors  in  the  determination  of 
diffraction  data, 

2.  errors  in  normalization  of  the  diffraction  data, 
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3. 


errors  in  the  tabulated  values  of  atomic 


scattering  factors  used  at  various  stages 
during  the  calculation  of  the  experimental 
PFD,  and 

4.  errors  in  the  tabulated  values  of  Compton 
scattering . 

Let  the  required  correction  in  the  function  F(k)  due 
to  the  above-mentioned  errors  be  AF(k).  The  Kaplow  correc- 
tion procedure  obtains  the  AF(k)  for  a given  set  of  experi- 
mental intensity  values  I (k) . This  is  achieved  as  fol- 

1 exp 

lows. 

Equation  (A-2.2)  is  first  determined  as  a function  of 
r using  the  IeXp(k)  values.  This  will,  in  general,  result 
in  a G(r)  function  containing  ripples  in  the  region  of  r 
values  between  zero  and  r - , where  r . is  the  position  of 

the  first  peak  of  the  PFD  (the  position  closely  corresponds 
to  the  first  peak  of  the  G(r)  function  as  well).  By  defini- 
tion, G(r)  should  be  a straight  line  given  by  -4TTrpo  in 
this  region.  The  error  in  the  G(r)  function  in  this  region 
is  the  difference  between  the  experimentally  determined  valu 
cf  G(r)  and  the  straight-line  value  in  this  region  and  is  re 

ferred  to  as  AG.  AF(k)  is  obtained  by  transforming  AG  and 

2 

multiplying  the  transform  by  I (k)/  Z f-  [Kaplow,  Strong 

P u.c. 

and  Averbach  (1965)].  The  true  value  of  F(k)  is  then  ob- 
tained by  subtracting  the  AF(k)  values  from  the  correspond- 
ing experimental  F(k)  values.  From  the  corrected  values  of 
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F(k)  the  corrected  intensity  values  can  be  obtained  as  a 
function  of  k using  Eq.  (A- 2.1).  These  corrected  values 
are  then  used  to  calculate  the  i(k)  of  the  right-hand  side 
of  Eq.  (3.2).  It  must  be  noted  that  AG  also  contains  a 
small  contribution  from  termination  ripples  but  the  effect 
of  these  ripples  had  to  be  ignored  in  our  present  investi- 
gation. An  attempt  was  made  to  eliminate  the  effect  of 
these  ripples  by  going  through  a cyclic  process,  but  the 
termination  ripple  values  were  too  small  and  hence  introduced 
error  ripples  in  the  correction  procedure.  The  final  PFD  is 
the  result  of  a Fourier  transform,  and  hence  contains  termi- 
nation ripples  due  to  the  termination  of  the  intensity 
values  at  a finite  S value.  However,  these  ripples  do  not 
cause  any  harm  because  they  are  also  present  in  the  calcu- 
lated pair  functions  (Appendix  3)  and  hence  in  the  calculated 
PFD. 

In  practice,  the  average  atomic  density  distribution 
is  expressed  in  terms  of  a RHO  given  by 


RHO 


o 


( ^ 
u.  c . 


I 

u.  c. 


V 


2 


Each  of  the  functions,  F(k),  AG  and  k"ti(k)e  a ^ are  sub- 
jected to  Fourier  transform  operations  at  various  stages 
of  the  Kaplow  correction  procedure.  These  operations  are 


carried  out  by  a computer  program.  An  important  feature 
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of  the  program  is  a subroutine  called  TERM.  The  subroutine 


TERM  is  designed  to  accomplish  the  following. 

2. 


■a2k2 


F(k),  G(r)  and  k i(k)e  are  functions  which  oscil- 

late about  the  value  of  zero.  The  last  available  values  of 
these  functions  may  or  may  not  be  zero.  If  the  last  value 
of  a function  is  finite,  unwanted  ripples  are  introduced 
during  the  transformation  of  the  function  due  to  this  factor 

alone.  The  subroutine  TERM  starts  with  the  last  value  of 

? -a2k“ 

the  functions  F(k),  AG  and  k~i(k)e  and  compares  it  to 

zero.  If  the  value  is  not  zero  it  takes  points  going  down 
from  the  last  value  until  the  function  crosses  the  zero 
value.  Then  the  subroutine  TERM  terminates  the  function  at 
that  point  after  assigning  the  value  of  zero  to  the  function 
at  that  point.  The  value  of  k at  which  the  function 
kZi(k)e  a ‘ is  terminated  (km)  is  also  given  in  the  print- 
out. The  same  value  of  k is  used  in  the  calculation  of 

m 

the  pair  functions  (Appendix  3)  used  to  calculate  the  theo- 
retical PFD. 


APPENDIX  3 


PAIR  FUNCTIONS 


The  pair  functions  are  defined  by  the  equation 


o 


sin  r.  .k  sin  rk  dk 
ij 


(A- 3.1) 


where  "i”  and  ”j”  stand  for  either  a Si  or  an  0 atom,  k is 

the  maximum  value  of  k used  in  the  calculation  of  the  experi 

mental  PFD.  The  value  of  k is  determined  by  the  subroutine 

m 

TERM  used  in  the  Kaplow  correction  procedure  (see  Appendix 
2).  The  value  of  a is  the  same  as  used  in  the  corresponding 
PFD  calculation.  The  other  symbols  have  been  already  de- 
fined in  Chapter  III. 

In  practice  it  is  necessary  to  only  calculate  the  auxil 
iary  pair  functions  given  by 


q:  . (x)  = 1/2  / ■ , * e'a  k cos xk  dk 

13  ° g (k) 


(A-3.2) 


and  the  pair  functions  are  given  by  the  relation 


(A-3.3) 
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Q^.fr+r^)  is  usually  negligible.  The  integration  involved 
in  Eq.  (A- 3. 2)  is  performed  by  a computer  program  which s in 
effect,  performs  a cosine  transform  of  the  function 


f.f . 


gUk) 


2,  2 

-a  k 


over  the  range  0 to  k^.  The  program  gives  as  output  a plot 
of  the  auxiliary  pair  functions  as  a function  of  x from  0 to 
2 A,  as  well  as  a punched  card  output  to  be  fed  into  the 
program  which  carries  out  the  calculation  of  the  theoreti- 
cal PFD.  A typical  set  of  auxiliary  pair  functions  used  in 
the  present  investigations  is  shown  in  Fig.  A-3.1. 


130 


Figure  A-3.1. 


The  auxiliary  pair  functions  Qij (x)  for 
the  three  combinations  Si-Si,  Si-O,  and  O-O. 


APPENDIX  4 


SILICIC  ACID 

A sample  of  silicic  acid  containing  88  percent  SiC^ 
and  the  remainder  H-,0  was  also  studied  during  the  course  of 
the  present:  investigation.  This  silicic  acid  can  be  repre- 
sented by  the  formula  SiC^'-j-y  ^0.  The  presence  of  half  a 
mole  of  water  per  each  SiC^  molecule  might  be  expected  to 
have  a pronounced  effect  on  the  structure.  Hence  a PFD 
study  was  made  of  this  silicic  acid. 

The  silicic  acid  sample  in  powder  form  was  packed  in 
the  pow'der  diffraction  jig,  and  the  Cu  and  Mo  diffraction 
data  obtained  using  methods  already  described  (Chapter  IV) . 
The  Cu  diffraction  data,  corrected  for  air-scattering,  back- 
ground noise  and  polarization  factor , were  multiplied  by  a 
constant  so  as  to  oscillate  equally  about  the  independent 
scattering  curve  (modified  and  unmodified  scattering  per 
unit  of  composition)  at  S values  near  1.2  X . The  Compton 
scattering  determined  from  the  International  Tables  for  X-ray 
Crys  tallography , Vol . Ill  (1962),  was  subtracted  from  the 
normalized  Cu  data  to  obtain  the  unmodified  scattering  up 
to  S values  of  1.2  X"1.  The  low  angle  data  were  linearly 
extrapolated  to  zero  at  S = 0 (shown  as  dashed  line  in  the 
unmodified  scattering  curve). 
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The  Mo  data  were  corrected  for  air- scattering , beam- 
off  background,  polarization  and  multiple  scattering  and 
normalized  assuming  complete  elimination  of  Compton  scatter- 
ing at  S values  greater  than  1.2  A . The  normalized  scat- 
tered intensity  in  absolute  units  per  unit  of  composition, 
Si02  — H90,  of  the  sample  of  silicic  acid  used  in  the  pres- 
ent investigation  is  shown  in  Fig.  A-4.1.  The  independent 
unmodified  scattering  intensity,  as  well  as  the  modified 
scattering  (up  to  S equal  to  1.2  A"1), is  also  shown.  Except 
for  a broadening  of  the  first  diffraction  peak,  the  intensity 
curve  for  the  silicic  acid  under  consideration  is  remarkably 
similar  to  the  intensity  curve  for  vitreous  silica.  The 
values  of  the  intensities  are  slightly  higher  for  the  case 
of  silicic  acid  but  this  can  be  explained  as  due  to  the 
greater  number  of  atoms  in  the  unit  of  composition  of 
silicic  acid. 

The  experimental  PFD  of  silicic  acid  is  shown  compared 
to  the  experimental  PFD  of  "T-D"  clear  fused  quartz  in  Fig. 
A-4.2.  Pronounced  differences  are  seen  in  the  case  of  the 
first  three  peaks  of  the  two  PFD.  The  first  peak  of  the 
silicic  acid  PFD  is  much  smaller  and  spread  out  than  the 
first  peak  of  "T-D"  clear  fused  quartz.  The  second  and 
third  peaks  seem  to  have  merged  together  with  a pronounced 
decrease  in  the  magnitude  of  the  third  peak.  The  fourth 
and  subsequent  peaks  of  the  silicic  acid  PFD  are  much 
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Figure  A- 4. 2.  The  experimental  PFD  for  "T-D"  clear  fused  quartz  and  silicic  acid. 
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smaller  in  value  than  the  corresponding  peaks  of  "T-D" 
clear  fused  quartz  but  their  distribution  seems  to  be  the 
same , 

It  is  known  that  water  in  silicic  acid  exists  in  com- 
bined form  as  OH  groups  rupturing  neighboring  Si -0-Si  bonds.” 
The  replacement  of  oxygen  by  a hydroxyl  group  may  be  ex- 
pected to  increase  the  Si-0  distances  involved  and  thus 
would  directly  affect  the  first  peak  of  the  DDF.  Also,  it 
would  be  expected  that  the  Si-lst  Si  interaction  in  silicic 
acid  should  be  much  weaker  than  that  in  vitreous  silica. 

This  explains  the  reduction  in  the  Si-1'  Si  contribution 
to  the  experimental  PFD.  However,  the  distribution  of  atoms 
starting  with  the  Si -2nd  0 onwards  should  not  be  affected. 
This  seems  to  be  the  case. 


*Sosman  (1965)  has  given  a reaction  to  explain  the  forma- 
tion of  silicic  acid  when  fine-ground  silica  W is  suspended 
in  CO2  free  water.  This  reaction  can  be  modified  for  the 
case  of  the  silicic  acid  under  consideration  in  this  study 
as  follows. 


°X  A A /° 

Si  Si  Si 

0/  V V \> 


+ J H2° 


OH  0 OH  HO  0 

V'  V NSi 

0'  V x 


0 


It  would  appear  that  only  three  out  of  every  four  Si-O-Si 
bonds  would  remain  intact  after  the  incorporation  of  water 
in  the  structure.  This  would  explain  the  decrease  in  the 
third  peak  of  the  silicic  acid  PFD. 
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